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Abstract 


Observations of 1 to 20 MeV/nuc oxygen, sodium, and siUfur ions in 
the Jovian magnetosphere are reported. Measurements made by the 
Cosmic Ray Subsystem on Voyager 1 and 2 are used to calcxilate 
abundances and energy spectra in the region from 5 to 20 Jovian radii (Rj). 
The phase space density of the oxygen ions calculated from the spectra 
has a positive radial gradient between 6 and 17 Rj, indicating an inward 
diffusive flow. The diffusion coefficient upper limit at 9 Rj is -10”® s”*. 
This limit, combined with the analysis of Voyager plasma observations by 
Siscoe cf qL. [198lj, implies an upper limit to the mass loading rate near lo 
of -10*® ions/s. 

The energetic oxygen lifetime is within an order of magnitude of the 
strong pitch-angle diffusion lifetime in this region, with the largest total 
number of particles lost between 7.5 and 12.5 Rj. It is shown that the 
losses are not due to geometric absorption by lo. absorption by dust 
grains, or energy loss in the plasma of the inner magnetosphere, and it is 
therefore postulated that the primary loss mechanism is pitch-angle 
scattering into the loss cone. 

The power delivered to the Jovian atmosphere by oxygen and sulfur 
ions with magnetic moments greater than 70 MeV/nuc-G (EiO.3 MeV/nuc at 
10 Rj), scattered into the loss cone between 6 and 17 Rj, is -5x10** W; a 
power comparable to the -1.2x10*® W required to produce the auroral 
emission observed by the ultraviolet spectrometer on Voyager [Broadfoot ef 
oL., 1981]. The measured dependence of the input power on magnetic 
moment threshold indicates that additional power is contributed by oxygen 
and sulfur ions with magnetic moments less than 70 MeV/nuc-G. Since 
most of the particles are lost between 7.5 and 12.5 Rj. the latitudinal zone 
in which the most energy is deposited in the atmosphere is —67® to -71®. 


V. 


Table of Conienta 


Acknowledgements 1 1 

Abstract i v 

1. Introduction 1 

2. The Experiment and Data Analysis 11 

2.1 Encounters with Jupiter 11 

2.2 The Instrument 16 

2.3 The Analysis Technique 22 

2.3.1 Element Identification and Energy Determination 22 

2.3.2 Intensity Determinations 28 

3. Observations 33 

3.1 Intensity Profiles 33 

3.2 Differential Energy Spectra 39 

3.3 Elemental Abundances 44 

3.4 Integral Energy Spectra 51 

3.5 The Phase Space Density 56 

3.5.1 Intensity - Phase Space Density Relationships 56 

3.5.2 Pitch Angle Distribution Corrections 58 

3.5.3 Phase Space Density Determination 63 

4. Discussion 73 

4.1 Overview 73 

4.2 The Oxygen Flow Pattern and Acceleration Process 74 

4.3 Diffusion Coefficient Determinations 75 

4.3.1 Lossless Diffusion 76 

4.3.2 Lossy Diffusion 80 



vi. 

4.4 Diffusive Flow Rates 84 

4.5 Particle Loss Mechanisms * 89 

4.5.1 Geometric Absorption by lo 89 

4.5.2 Absorption by Dust 90 

4.5.3 Oxygen Charge State 91 

4.5.4 Energy Loss in Plasma 92 

4.6 Energy Flow Rates 95 

5. Conclusions 104 

Appendix A - Definition of Terms 107 

Appendix B - High Counting Rate Calibrations 109 


References 


118 


1 


Chapter 1 
Iniroductioii 

Information provided by the Voyager spacecraft during their 
encounters with Jupiter in 1979 has changed our basic understanding of 
the particle composition and energetics of the Jovian magnetosphere. 
Oxygen and sulfur ions are now known to play important roles in all 
particle populations, from the plasma near lo to the MeV per nucleon 
(MeV/nuc) ions of the radiation belts. 

The source of these ions is the Galilean satellite lo. whose 
volcanoes eject SO 2 onto its surface and into its atmosphere. Some of 
the gas that escapes from lo becomes ionized and is picked up by the 
rapidly rotating magnetic field of Jupiter (9^ 55“ period). The energy 
produced by this acceleration process helps heat the torus of plasma 
that encircles Jupiter at the orbit of lo, producing intense ultraviolet 
radiation. As the plasma diffuses outward from the torus, it is pulled into 
a thin equatorial disk by the centrifugal forces acting on the heavy ions. 
Once the ions reach ths outer magnetosphere, a mechanism that is not 
yet fully understood further accelerates some of the particles, producing 
a source of energetic oxygen and sulfur ions. These ions then diffuse 
away from the region, both outward along open field lines into the Jovian 
magnetotail producing an oxygen and sulfur wind, and back inward 
towards lo. 

The inwardly diffusing ions gain energy as they move toward 
stronger magnetic fields, till they acquire energies in excess of an MeV 
per nucleon. Although they are small in number compared with the ions 
in the plasma, the total energy they bring into the inner magnetosphere 
is significant. This thesis concerns observations naade of the energetic 
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oxygen and sulfur ions inside 20 Jovian radii (Rj; 1 Rj- 7.14x10* km) by 
the cosmic ray instrument on Voyager 1 and 2. To put these 
observations in context, we will now briefly discuss the measurements 
made of the Jovian magnetosphere prior to the Voyager encounters, and 
describe in more detail the relevant results from Voyager. 

In the 1950s Jupiter was discovered to be a source of nonthermal 
radio emission at both decametric [Burke and Franklin. 1955] and 
decimetric [Drake and Hvatum, 1959; Berge and Gulkis. 1976] 
wavelengths. It was found that the decametric radiation is modulated by 
Io*s postion [Bigg. 1964], and postulated tliat the emission occurs in the 
ionosphere of Jupiter at the feet of the fleld lines that pass through the 
satellite. Equating the observed cutoff frequency of the emission with the 
local electron g 3 rrofrequency led to an estimate of 7 G for the equatorial 
surface fleld strength of Jupiter [Goldreich and Lynden-Bell, 1969]. The 
decimetric radiation on the other hand was interpreted as synchrotron 
emission from electrons trapped in a dipolar magnetic fleld [Field. 1959]. 
Hence, even before any in situ measurements had been made of the 
Jovian magnetosphere, it was known that the planet had a strong 
magnetic fleld, containing trapped energetic electrons. 

The flrst spacecraft to make direct measurements o! the Jovian 
magnetosphere were Pioneer 10 and 1 1 , which flew by Jupiter in 1 973 and 
1974 respectively. The instruments on the spacecraft include a 
magnetometer (two on Pioneer 11). a plasma analyzer, an ultraviolet 
photometer, and four particle detectors covering the energy range from 
~100 keV to ~50 MeV for electrons and ions. The results from each 
instrument can be foimd in three collections of articles [Snencs, 183, 
301-324, 1974 (Pioneer 10); Science, 188, 445-477. 1975 (Pioneer 11); J. 
Gsophys. Rbs.. 79, 3467-3694. 1974], with interpretations of these results 
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and reviews of Jovian research through 1975 contained in the book 
«Aiptfer [Gehrels, 1976]. In addition, there are four review articles on the 
Jovian magnetosphere baaed on Pioneer observations [Coroniti, 1975; 
Kennel and Coroniti. 1977; Schulz. 1979; Goertz and Thomsen. 1979]. 

The Pioneer results that pertain to the present discussion can be 
summarized as follows: The magnetic field is basically dipolar, with on 
equatorial strength of 4 G and an axis tilted 1 1^ with respect to Jupiter’s 
spin axis. However, in the region outside ~15 I^. and particularly outside 
30 Rj, the field was found to be distorted in the vicinity of the equatorial 
plane by an additional radial component. The distortion was postulated 
to be due to a warped current sheet caused by centrifugal forces acting 
on low energy magnetospheric plasma [Smith ef ol.. 1974; Van Allen et al., 
1974]. 

High intensities of energetic electrons and protons are trapped in 
the Jovian magnetosphere. Although the measured intensities above a 
given energy threshold generally increase toward smaller radial distances, 
the phase space densities of the electrons and protons at constant 
particle magnetic moment. M, (M«Pj^^/B) were found to iccrease inward 
[McQwain and Fillius, 1975; Baker and Goertz. 1976; Thomsen at ai., 1977]. 
The phase space density at constant M is the quantity that controls the 
magnitude and direction of particle diffusion (see e.g. Schulz and 
Lanzerotti [1974]). so that the observations imply an inward radial 
diffusion of energetic electrons and protons. Particles conserving their 
magnetic moments gain energy as they diffuse toward stronger magnetic 
fields. In the Earth's magnetosphere, for example, this process is the 
dominant mechanism for producing radiation belt particles [Schulz and 
Lanzerotti, 1974. and references therein]. However, in the Jovian 
magnetosphere the combination of large numbers of ibIieV electrons in 
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the outer magnetosphere [Gocrtz, 1978] and >50 MeV electrons (M = 2xl0* 
^ MeV/G at 3 I^) in the inner magnetosphere [FUlius and McUwain. 197< 

r 

Coronitl, 1975] requires an additional acceleration mechanism [Sentman 
€< at., 1975; Carbary 1976; Goertz, 1978 with numerical modelling by 
Borovsky et at.. 1981]. 

Loss mechanisms act on the energetic particles as they diffuse 
inward. Sharp decreases in electron and proton counting rates at the 
orbit ol Europe. lo. and Amalthea and changes in the pitch-angle 
distributions across the orbit of lo indicated that some of the particles 
are absorbed by the Jovian moons [Thomsen, 1979], There is also strong 
evidence that additional losses occur in the inner magnetosphere 
[Barbosa and Coroniti, 1976; Goertz et at., 1979], most likely caused by 
pitch- angle scattering into the loss cone [Fillius ct at., 1976]. 

The only direct indication of the presence of heavy ions in the 
magnetosphere from the Pioneer encounters was from the University of 
Chicago fission cell detector, designed primarily to measure proton 
intensities but also sensitive to heavier nuclei. Signal characteristics in 
the inner magnetosphere were inconsistent with a proton dominated 
response, leading Simpson ef at. [1974, 1975] to suggest that nuclei with 
nuclear charge, Z, greater than one were contributing to the counting 
rate. It has recently been shown using Voyager measurements that it is 
possible the fission cell was counting mostly oxygen and sulfur ions in the 
inner magnetosphere [Gehrels at at... 1961]. Also, it is now known that 
the plasma analyzer on Pioneer 10 was detecting oxygen and sulfur ions 
near lo [Intriligator and Miller, 1981], although data from the instrument 
were initially analyzed assuming a proton plasma [Frank at at., 1976]. 

Conclusive evidence of heav)' ions and neutrals in the 
magnetosphere was provided prior to the Voyager encounters by ground- 
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baaed optical observations. In 1972. Brown [1974] observed emission 
Irom a cloud of neutral sodium atoms surrounding lo. and it was 
predicted that sodium plasma and energetic ions would be present in the 
magnetosphere [Eviatar tt al., 1976]. Emissions from singly ionized sulfur 
were discovered in 1975 [Kupo et al., 1976], and since that time many 
other neutral and ionized species have been observed, including KI. 01, 
on, Oin, Sni, and SIV (see review by Pilcher and Strobel [1981]). 
Although some consideration was given to the effect heavy ions would 
have on the magnetosphere [HiU and Michel. 1976], the significant role 
they play at all energies was not appreciated until the Voyager 
encounters. 

Voyager 1 and 2 encountered the Jovian system in March and July 
of 1979. The instruments on the spacecraft include a magnetometer, a 
plasma analyzer, a plasma wave detector (10 Hz * 56 kHz), a planetary 
radio experin^nt (1.2 kHz - 41 MHz), an ultraviolet spectrometer, a low 
energy particle experiment (electrons 14 keV •• 10 MeV, ions 30 keV - 150 
MeV with elemental identification for 2200 keV), and a cosmic ray 
experiment (electrons 3—100 MeV, protons -iron >^l-«-20 MeV/nuc 
depending on element). There are several collections of articles 
concerning Voyager ren'ilts and interpretations [5bwnc«, 2M. 945-1008. 
1979 Cv^oyager 1); Science, 206. 925-996, 1979 (Voyager 2); Nature, 280, 
725-806, 1979; Ceophys. Res. Lett., 7, 1-66, 1980; J. Geophys. Res., 86, 
1981]. The discussion of the results relevant to this thesis will be divided 
into the following topics: energetic particles, the lo plasma torus, the 
aurora, and the particle and energy budget 

Energetic particles - In the outer magnetosphere, the 
composition of ions with Z>2 between 7 and 14 MeV/nuc is similar to that 
of solar energetic particles [Vogt et al., 1979a, bj. However, inside -12 Rj, 
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the coinposiUon is dominated by oxygen and sul!ur. In the inner region 
the sulfur to oxygen ratio at constant energy per nucleon is between 0.6 
and 0.8. A similar change in composition occurs for ions with energies in 
Uie range 0.60 to i.l5 MeV/nuc. except that, in the outer magnetosphere, 
sulfur is found to be approximately ten times more abundant relative to 
helium than in solar flare events [Krimigis tt al„ 1979; Hamilton at al., 
1961]. As in the case of energetic electrons and protons observed by 
Pioneer, the phase space densities of the energetic oxygen ions [Vogt at 
al., 1979b; Gehrels at al., 1961] and protons [Armstrong at al., 1981] 
detected by Voyr^er decrease toward smaller radial distances indicating 
inward diffusion. j>eyond ~150 Rj in the magnetotail, ions with energies 
of -100 keV stream outward in the anti-solar direction. This 
magnetospherio wind is rich in oxygen and sulfur, indicating that lo is 
the source of the particles. 

lo plasma torus - The two most important Voyager discoveries 
concerning the magnetosphere were the intense Jovian aurora and the 
torus of hot plasma composed predominantly of oxygen and sulfur ions 
that encircles Jupiter at the orbit of lo. The source of the ions in the 
torus is lo's SO 2 atmosphere fed by volcanic activity on the satellite. The 
torus is divided into two parts, a cold (ion terr^ieratures <1 cV) inner 
region inside 5.5 R; with a scale height of -0.? Rj from the centrifugal 
equator, and a hot (ion temperature -30 eV) outer region between 5.7 
and -7 Ri with a scale height of -I Rj [Bagenal and Sullivan. 1981]. (The 
orbital position of lo is 5.9 R;.) The grotmd-based observation of SH Unes 
[Kupo at al., 1976] discussed earlier was of emission from the cold torus. 
The hot torus emits strongly in the ultraviolet, and was first observed by 
the ultraviolet instrument on Voyager 1 [Broadfoot at al., 1979]. It has 
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subsequently been studied in the ultraviolet by Voyager 2 [Sandel ef al., 
1979] and by the International Ultraviolet Explorer [Moos and Clarke, 
1981]. The total rate of radiation energy-loss from the hot torus is 
estimated to be 3x10^^ V [Broadfoot at al., 1979; Shemansky, 1980a]. 

The charge density in the torus due to species with A/Z*28. where 
A is the atomic mass number and Z* is the ion charge state, (i.e., due to 
oxygen and sulfur in low ionization states) was directly measured by the 
plasma instrument [Bagenal and Sullivan. 1981], and was foimd to be 
equal within statistics to the electron density inferred from 
measurements of the electron plasma and upper hybrid frequencies by 
the planetary radio experiment [Warwick at oL., 1979]. The implication is 
that ions with A/Z*<8 (Le. hydrogen) are a negligible component of the 
plasma. The charge density is >10® cm"® between 5.7 and *^7 Rj and has 
several peaks in this region (see Figure 4.3a of this thesis). The density 
remains high (i^lO® cm"®) out to ~9 Rj where the plasma disk begins to 
form. Plasma wave activity is enhanced inside ~20 Rj and is particularly 
intense in the region of the hot torus [Scarf at al., 1979; Gurnett at al., 
1979; Scarf at al., 1981]. 

Jovian aurora — Auroral emission (H Ly^ and Hg Lyman and 
Werner bands) was observed in the Jovian atmosphere at high latitudes 
near both poles by the ultraviolet spectrometer [Broadfoot at al., 1979; 
Sandel at al., 1979]. The total power required to produce the. observed 
luminosity is estimated to be 1.2x10'® W [Broadfoot at of., 1981] and is 
presumably supplied by the precipitation of magnetospheric particles. 
The auroral zone is approximately 6000 km wide from north to south 
corresponding to in latitude, and is located between 60® and 80° 
latitude [Broadfoot at al., 1981]. The location of the equatorward 
boundary of the zone is generally consistent with the magnetic projection 
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of the plasma torus onto the atmosphere [Broadfoot ef al., 1961], leading 
to the suggestion that electron precipitation from the torus drives the 
aurora [Sandel at al., 1979; Thome and Tsurutani, 1979]. There are, 
however, large differences in the positions of different aurorae and 
uncertainties in the positions of all of them due to the low polar viewing 
angle of the measurements, so that the magnetospheric location of the 
particle source is not constrained to the torus. Observations of optical 
auroral emission by the imaging instrument on Voyager [Smith at al., 
1979: Cook and Jones. 1961] may be useful in mapping the location and 
depth of the auroral activity, thereby indicating the energy and t)rpe of 
particle that is causing the emission at various depths and latitudes. 

Particle and energy budget - There has been wide interest 
[Thorne and Tsurutani, 1979; Goertz, 1980; Eviatar and Siscoe, 1980; 
Shemansky, 1900b; Thome, 1981; Sullivan and Siscoe, 1901] in determing 
the source of power that heats the torus (~3xl0** W) and produces the 
auroral activity (~1.2xl0^^ W). A possible energy source is the 
acceleration of newly ionized oxygen and sulfur ions near lo to corotation 
speeds by Jupiter's magnetic held. The ions acquire a longitudinal 
guiding-center velocity equal to the corotation velocity of the magnetic 
field at Id’s orbit, and. assuming no local retardation of the corotation 
velocity near lo. a gyrospeed equal to the difference between lo's orbit 
velocity and the corotation velocity. This difference is 66 km/s, so that 
oxygen and sulfur ions acquire 262 and 524 eV respectively in the 
corotating frame. Therefore, the rate at which ions must be created and 
accelerated (the mass loading rate) in order to heat the torus is *^5x10^ 
ions/s (assuming twc oxygen ions for every sulfur ion due to dissociation 
of SOg) and to produce the auroral power is '^2x10^ ions/s. However, 
Shemansky [1900b] has shown that excitation of ions during the 
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ionizaLion process would lead to ultxaviolei emissions of specific 
waveleniiths from the ionization region, and that no such emissions were 
detected in the vicinity of lo. The upper limit to the mass loading rate 
near lo consistent with no oberved radiation is <^10*^ ions/s. The 
measurement does not rule out the possibility that ionization occurs in 
other parts of the torus, away from lo. but there is another constraint 
that limits this case. 

The observed abundances of the different ionization states of 
sulfur in the hot torus during Voyager 1 encounter require that the ion 
residence time in the torus is on the order of 100 days [Shemansky. 
19fl0b]. The total number of ions in the hot torus is ~5xl0^ [Sullivan 
and Siscoe. 1961], so that this analysis gives a mass loading rate on the 
order of <^6x10^ ions/s. Although this estimate is not accurate to a 
factor of two or probably even four, it does suggest that the mass loading 
rate is no lower than ~10*’ ions/s, and no higher than a few times 10*® 
lons/s. The mass loading rate of 5x10*® ions/s required to heat the 
torus is therefore somewhat high, and the 2x10*® ions/s required to 
power the aurora is certainly unrealistic. This problem is often referred 
to as the energy crisis of the inner magnetosphere. 

A possible solution to the auroral power requirement is to 
precipitate inwardly diffusing energetic protons [Goertz, 1900] or heavy 
ions [Thome, 1981] instead of torus electrons. In this case the auroral 
energy is supplied by the process that accelerates the energetic particles. 
Thorne [1961] also points out that some of the secondary electrons 
created in the atmosphere by heavy ion precipitation could escape back 
into the magnetosphere and help heat the torus. One of the goals of the 
present analysis is to determine the amount of power delivered to the 
Jovian atmosphere by energetic oxygen and sulfur ions lost from the 
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inner magnetosphere. 

This thesis is based on a study of energetic Z>2 ions in the Jovian 
magnetosphere. In previous reports [Vogt mt of.. 1979a.b; Gehrels ef of., 
1981], the elemental composition of the energetic ions, the oxygen energy 
spectrum, and the phase space density profile between 6 and 17 Hj for 
680 MeV/nuc-G oxygen ions were determined. In this report oxygen, 
sulfur, and sodium energy spectra are obtained, and the phase space 
density profile is calculated with improved radial resolution for several 
oxygen magnetic moments. These profiles are used to specify limits on 
the diffusion coeflfcient in the inner magnetosphere. The inward diffusive 
flow rate of particles and energy is then determined and estimates made 
of the power delivered by the oxygen and sulfixr ions to the Jovian 
atmosphere. 



11 


Chapier 2 

The Experiment and Data Analysis 
2.1 Encounien vith Jupiter 

The data used in this study were obtained with the Low Energy 
Telescopes (LETs) ol the Cosmic Ray Subsystem (CRS) on board Voyager 1 
and 2 during their encounters with Jupiter in 1979. Trajectory and 
timing information for both encounters are given elsewhere [Sbience 204, 
945-1008. 1979; Sbisnce 206, 925-996. 1979] and will be only summarized 
here. Figures 2.1 and 2.2 show the spacecraft trajectories with respect to 
the Jovian magnetosphere and satellites. Typical positions of the bow 
shock and magnetopause boundaries (see Appendix A for definitions of 
roagnetospheric terms) that delimit the magnetosphere are indicated in 
Figure 2.1. These positions can. however, vary by as much as ±30 Rj on a 
time scale of hours. 

From February 20 through March 2. 1979, Voyager 1 crossed the 
bow shock and magnetopause a number of times as these bouiidaries 
convected back and forth past the spacecraft (exact crossing Uxnes are 
given by Lepping at al. [1980]). The final entry into the magnetosphere 
was at 0227 universal time (UT, spacecraft event time) on March 3 at a 
zenocentric distance of 47 I^. Closest approach to Jupiter was at 1205 
UT on March 5 (day 64 of 1979) at a distance of 4.9 Rj. just inside the 
orbit of lo (Figure 2.2a). With a periapsis near the orbit of lo, a 
maximum amount of time was spent observing the satellite and its 
interaction with the magnetosphere. The spacecraft flew almost directly 
underneath lo on the outbound crossing of the orbit, allowing the current 
system in the magnetic flux tube that connects it and Jupiter to be 
studied. Voyager 1 crossed the magnetopause outbound several times on 
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FS^rure 2.1 

Voyager 1 (VI) and 2 (V2) trajectories through the Jovian magnetosphere 
as viewed from the north. Jupiter is at the origin of the coordinate 
system. Typical bow shock and magnetopause positions are shown. 
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Figure 2J3 

Trajectories of a) Voyager 1 and b) Voyagei 2 with respect to the 
orbits of Jupiter’s moons, as viewed from the north. 
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March 15 at a distance of ~160 I^. 

The Voyager 2 bow shock and magnetopause crossings occurred 
between July 2 and 5« 1979. with final entiy into the magnetosphere at 
1640 UT on July 5 at a distance of 62 Rj. Closest approach was at 2229 
UT on July 9 (day 190) at a distance of 10.1 just outside the orbit of 
Europa (Figure 2.2b). The Voyager 2 outbound trajectory was the deepest 
penetration of the Jovian magnetotail of any spacecraft to date, with the 
final magnetopause crossing occurring at a distance of ~260 on August 
3. 

Both spacecraft had inbound and outbound trajectories that were 
within latitude of Jupiter's spin equatorial plane. Therefore, since 
the planet's magnetic dipole axis is tilted 10.6” with respect tc its spin 
axis, the spacecraft passed through the magnetic equatorial plane twice 
per planetary rotation period of 9^ 55® (expect beyond «-0O Rj in the 
magnetotail where the magnetic equatorial plane is bent toward the spin 
equatorial plane). This variation in spacecraft magnetic latitude caused a 
periodic modulation in the measured particle intensities due to the fact 
that the intensities of trapped energetic particles are organized about 
the magnetic equatorial plane, which is their bounce S 3 rmmeti 7 surface. 

2.2 The Instrument 

This section describes the features of the LETs that are important 
for the present anal 3 rsis. General descriptions of the CRS instrument and 
electronics are given by Stone al. [1977] and StilwcU al. [1879]. 
Detailed inlormaticn about the CRS instrument such as discriminator 
settings, data formats, and rate specifications are in the Voyager Science 
Requirements Document [Garrard, 1976]. 
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The LETT system provides measurements of particle intensities, 
energies, and elemental composition in the lieV/nucleon energy range. It 
consists of four nominally identical particle telescopes. LET A, B, C. and 
D. oriented to allow detection of anisotropies in the particle arrival 
directions. During the Voyager 1 and 2 Jovian encounters. LET A and C 
were turned off inside ~20 Rj to reduce the noise level at the input of the 
pulse-height analyzers shared with LET B and D respectively. Therefore, 
since the analysis that follows concentrates on the region inside 20 Rj. 
only data from LET B and D will be used. 

A schematic diagram of a LET is shown in E^ure 2.3. The telescope 
is a stack of four solid-state silicon surface-barrier detectors (LI. L2, L3. 
and L4). each of which is a thin circular disk of cm radius. A 3 /im 
thick aluminum window provides light shielding for LI. When adequate 
bias is applied between the thin ("-40 /ug/cm?) metal contacts on each 
surface of the detector, the silicon becomes fully depleted. 

A charged particle penetiating or stopping in the detector creates 
electron-hole pairs along its path that collect on the contacts, producing 
a small current pulse that is fed into a charge sensitive preamplifier and 
a shaping postemplifler. The height of the resulting output pulse is 
proportional to the particle energy loss in the detector. The signal from 
each detector aii^lifler is then compared with a threshold level in a 
discriminator circuit, and, if a specified coincidence requirement between 
the four detectors is met, the signals from LI, 12 , and L3 are pulse-height 
analyzed. Coincidence requirements are specified by sending commands 
to the instrunaent. 

The data used in the present study were obtained with Voyager I 
in the L1L2 command itate in which a particle is required to exceed the 
LI and L2 discriminator thresholds but no requirement is made on L3, 
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Figure 2.3 

Schematic diagram of a LET. The detector areas and thicknesses that 
are given are nominal values. Laboratory measurements of the flight 
detector thicknesses and active areas [Gehrels emd Cummings^ 1960] were 
used in the analyst. 
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and with Voyager 2 in the L1*L2 L3 command state in which a particle is 
required to exceed the discriminator thresholds in LI. L2, and L3. During 
the inteiplanetary cruise parts of the missicn. L4 is used in 
anticoincidence in these two command states to eliminate penetrating 
particles. However. L4 was turned off inside ~20 Rj in the Jovian 
magnetosphere to minimize accidental event rejection due to vhe high 
radiation environment. The LET geometry factors for pulse-height 
analyzed events (PHA events), and the energy ranges over which oxygen 
and sulfur nuclei can be measured are listed in Table 2.1 for both 
command state;i. 

In addition to pulse heights, the instnunent also measures the 
rates at which discriminator threshold levels are exceeded and 
coincidence conditions are met The following rates are used in thi.s 
study: 

1) detector singles rates. The counting rates of the LI. L2, and L3 
discriminators; called Ki, Rs, and Rg. 

2) Z>2 singles rate. The counting rate of the Z>2 discriminator. The Z>2 
discriminator is triggered when a particle has energy losses El. £2, and 
E3 in LI. L2, and L3 such that El 4 - 0.42 £2 4 - 0.20 E3 > 0.6MeV. It is called 
the Z>2 discriminator because neither protons nor alpha particles within 
the acceptance cone of the telescope can satisfy this condition. Note 
that the Z>2 rate does not require energy loss in all three detectors (i.e. 
it is not a coincidence rate), and therefore, since the geometry factor 
defined by the collimator and LI (4.6 cm^ sr) is much larger than that 
defined by the collimator and L2 or L3 (0.5 cn^ sr). it is primarily the 
rate of Z>2 nuclei losing more than 9.6 UeV in LI. The geometry factor 
and oxygen and sulfur energy ranges for Ute Z>2 rate, aie listed in Table 
^l. 
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Table 2.1 

Nominal LET Parameters 



Command 

State 

Geometry * 
Factor 

2 

(cm -sr) 

Element 

Energy + 
Range 

(MeV/nuc) 

event analysis 

L1-L2 

0.44 

oxygen 

sulfur 

3.6 - ^-23. 
4.5 - 'vSO, 


L1-L2-L3 

0.44 

oxygen 

sulfur 

6.0 - 'V.23. 
7.5 - -vSO. 

Z>2 rate 

- 

4.6 

oxygen 

sulfur 

'ul.l - 25. § 

-x-C.g - 120. 


* Each of four telescopes. 

T Energy loss in the 3 ym aluminuin window has been added to 
the detector and Z>2 thresholds to give the incident energy 
of the particle. 

If Event analysis includes events that penetrate L3 (L4 off). 

The maximum energies are those at which the track foldback 
on an E2 vs E3 plot crosses the oxygen track for sulfur and 
the carbon track for oxygen. 

5 The dependence of the low energy threshold of the Z>2 rate on 

rate and pulse height is given in Table B.2. The oxygen thresholds 
used in the analysis are listed in Table 3.2, 
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3) the Z>2 coincidence rate. The counting rate of particles exceeding 
the discriminator thresholds in Li and L2 (L1^L2 command state) or in LI, 
L2, and L3 (Ll'L2*L3 comnw^d state), and in addition exceeding the Z>2 
threshold; called th^ LZ3 rate. The conditions for this rate are identical 
to those for pulse-height analysis. During the Jupiter encounters, the 
spacecraft data rate limited the rate at which Z>2 PHA events were 
received to approximately one per second (per LET). Since this was a 
small fraction of the incident rate of these events, the LZ3 rate was used 
to renormalize the intensities. 

2.3 The Analysis Technique 

The technique for converting rate and pulse height informaUvnn 
into elemental abundances, energy spectra, and absolute intensities will 
now be discussed. 

2.3.1 EHement Identification and Energy Determination 

A particle entering a LET is stripped of most of its electrons in the 
3 /Am aluminum window, so that no information is available about its 
charge state prior to entering the telescope. The LETs employ a AE — E 
technique to identify the nuclear charge of particles that penetrate LI 
and stop in L2 or L3. This technique and its application to the LETs is 
discussed in detail by Cook [i960], and will be only briefly summarized 
here. 

The energy loss, AE, of a nucleus with charge Z and velocity v in LI 
or L2 is roughly proportional to ZVv*. Therefore, by measuring a 
particle's energy loss in LI or L2 and its residual kinetic energy, E, as it 
stops in subsequent detectors, its cnarge can be determined. Figure 2.4a 
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ngure 2.4 

Energy-loss plots of Z>2 PHA events. The energy losses were calculated 
from the pulse heights returned by the instrument using laboratory 
calibration data [Povlis, 1900] and gain correction factors obtained from 
flight data [Cook, i960]. 

a) E2 vs E3 for all Voyager I LET D events from day 64/0936 to 64/1405 
UT (the 4.9 and 5.3 Rj regions in Table 3.1). The oxygen, sodium, and 
sulfur tracks are indicated. Events in the bottom left-hand corner have 
not been plotted. 

b) E2 vs E3 for all Voyager 2 LET B events from 190/1200 to 191/0800 
UT. Events in the bottom left-hand comer have not been plotted. 

c) El vs E2+E3 for the same events shown in figure a. Events along the 
El axis have not been plotted. 
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shows the L2 and L3 energy losses of Z>2 PHA events detected by LET D 
near lo. In this , case, E2 is the A£ measurement and E3 is the residual 
energy. Particles stopping in L3 are seen to form clear tracks 
corresponding to oxygen, sodium, and sulfur nucleL The events along the 
left-hand edge of the figure are due primarily to oxygen and sulfur 
events stopping in L2 that were in accidental coincidence with a Z=1 
particle (a proton or several electrons) that exceeded the L3 
discriminator threshold. The charge of these events can be determined 
using their LI and L2 pulse heights. There are a substantiad number of 
events in the bottom left-hand corner of the figure that have not been 
plotted. These comer events are due primarily to oxygen and sulfur 
nuclei stopping in or penetrating LI (exceeding the Z>2 threshold), but 
not hitting L2 or L3. The small pulse heights in L2 and L3 are due to Z=1 
particles in accidental coincidence with the heavy event. Figure 2.4b is a 
similar plot obtained near Voyager 2 closest approach (~10 Rj) by LET B. 

The width of the oxygen track in Figure 2.4a is a, = 0.2 charge 
units, compeared with the LET charge resolution of o, = 0.1 charge units 
(for Z near 8) measuied during solar energetic particle events in the 
interplanetary cruise part of the mission [Cook, 1980]. The poorer 
charge resolution obtained in the Jovian magnetosphere is due to 
eneigy-loss broadening caused by Z=1 particles in coincidence with the 
heavy ion. 

If the heavy ion rate itself becomes large enough in a detector that 
a significant number of the events are in comeidence with each other, a 
much larger broadening occurs. In L2 and L3 the spectrum and geometry 
factor combine to keep the oxygen and sulfur counting raie small 
(^10® s”'), so that this effect is not a problem in these detectors. 
However, as discussed in Appendix B, the oxygen and sulfur rate in LI 
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(b&sically the Z>2 rate) is high enough in some regions to cause 
substantial broadening; for instance, near lo ~30% of the oxygen and 
sulfur events are in coincidence with another heavy ion and are 
consequently displaced from their respective tracks. The displaced 
events form a backgroimd of events above the oxygen and sulfur tracks. 
This background can be clearly seen in Figvire 2.4c, which shows the 
energy losses El and E2+E3 for the same data as that shown in Figure 
2<4a. Fortunately, this problem can be avoided for events stopping in L3 
since element identification can be made using only E2 and E3. Almost 
all analysis in this study is based on such three-detector events. Only 
near Ic was the analysis extended to lower energies using the events 
stopping in L2 (two-detector events). For this case a pileup correction 
factor was applied to the two-detector portion of the spectrum, as 
described in Appendix B. 

For the three-detector events, the particle energy prior to entering 
the telescope was calculated from E2 and E3 using laboratory 
measurements of the LI thickness [Gehrels and Cummings, 1960] and the 
cn^irical range-energy relationship that Cook [i960] developed for the 
LETS. The telescope window thickness was assumed to be 3 ;um as 
specified by the nxinufacturer. For the two-detector events, the total 
energy is the sum of El and E2, plus a small term due to energy loss in 
the window. 

Intensity Detenninations 

There are two sources of information from the LETS concerning the 
ab^nlute Intensities of Z>2 particles: the Z>Z rate and the combination of 
PHA events and the LZ3 rate. As discussed in the instrument section, the 
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Z>2 rate is primarily the counting rate of nuclei with nuclear charge 
greater than two. losing more than 8.6 MeV in LI. Therefore, if the 
fraction of rate counts due to a certain element can be determined, the 
Z>2 rate can be used as an integral intensity measurement of that 
element. (The rate is not strictly speaking an integral measurement 
since particles above a certain energy are not counted, but spectra of 
interest in this study are all soft enough in the vicinity of the upper 
energy limit that the rate is effectively an integral.) 

The problem of determining this fraction is made easier in the 
inner Jovian r'agnetosphere by the fact that, in the MeV/nuc energy 
range, oxygen and sulfur are much more abundant than other species. 
Near lo, the piilse-height histogram of Li. shown in Figure 2.ba, reveals 
two separate peaks due to penetrating oxygen and sulfur nuclei. The 
events shown in this figure all have small 12 pulse heights caused, in 
roost cases, by a Z=1 particle in accidental coincidence with the heavy 
ion. The histogram therefore represents a sample of the Z>2 rate and 
can be used to directly measure the oxygen and sulfur fractions in the 
rate. In other regions of the magnetosphere, pulse-height histograms of 
Li did not have separate peaks as can be seen in the example from LET B 
near Voyager 2 closest approach shown in Figure 2.5b. In these regions 
the oxygen and sulfur fractions of the Z>2 rate were estimated from 
higher energy spectral data. 

The effective threshold level of the Z>2 discriminator can increase 
in high radiation environments due to baseline shifts at the outputs of 
the detector amplifiers. As described in Appendix B, this effect has been 
measured in the laboratory and is taken into account in all analysis. 
Also measured in the laboratory and described in the appendix is the 
hvedme correction for the Z>2 rate. 



yigure 2.5 


Histograms of the Li energy loss for Z>2 PHA events. The time periods 
and telescopes are the same as for a) Figure 2.4a and b) Figure 2.4b. 
Only events with E2 and E3 less than 15 MeV have been plotted so that 
these events correspond to those in the suppressed bottom left-hand 
corners in the previous figures. The triangles on the El axis indicate the 
Z>2 discriminator threshold level for LI during the time periods shown. 
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The PHA events are in general easier to convert into intensities 
than the Z>2 rate since the nuclear charge is measured for each event 
Also, the energy of each event is measured, so that differential as well as 
integral intensities can be determined. As discussed in the instrument 
section (2.2). the LZ3 rate is used to convert the PHA event rate into 
intensities. The formula for calculating the integral intensity of species k 
with energy greater than £' is 


Jk(E>E') = 


nph«-k(E>E0 

npk,^t 0*44 cm?-sr 


• LZ3 


( 2 . 1 ) 


where npk«-k(E>E') is the number of PHA events of species k with energy 
greater than E' measured in a given time interved. npk«-tot is the total 
number of PHA events in that time interval, and 0.44 cm?-sr is the 
geometry factor for event anal 3 rsis. The LZ3 rate used in equation (2.1) 
includes a livetime correction factor that was determined using 
laboratory measurements (Appendix B). 
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Chapter 3 
Obaenrationa 

Energetic particle observations made by the CRS instruments in 
the Jovian magnetosphere are summarized by Vogt at at. [1979a,b] and 
Schardt at al. [1961]. This study deals primarily with observations of Z>2 
ions inside ~20 Rj. The emphasis will be on oxygen ions, since oxygen 
was the only species observed with adequate statistics to allow a detailed 
investigation of the radial dependence of its intensity and energy 
spectrum. The abundances and spectra of carbon, sodium, and sulfur 
will, however, also be presented v lere statistics allow. 

3.1 Inteusitj Profiles 

Figure 3.1 shows the ireasurcd Z>2 intensities lor Voyagers 1 and 2 
inside ~20 Rj. The Z>2 intensity is determined from the livetime 
con^cted Z>2 rate and its geometry factor (4.6 cm? sr), and therefore, 
inside 20 Rj, is primarily the intensity of oxygen ions with energies 
greater than ~1.1 MeV/nuc and sulfur ions with energies greater than 
~0.9 MeV/nuc (Table 2.1). 

There is a pronounced decrease in the Z>2 intensity inside ~7 Rj 
on Voyager 1 due to particle losses in the vicinity of lo. The minima that 
occur in the intensities at 5.6 Rj inbound (day 64 0934 UT) and 5.4 Rj 
outbound (day 64 1406 UT) indicate an abrupt inner edge of the lossy 
region. Their location can therefore be compared with that of lo and 
that of plasma and electromagnetic wave features in the region to 
determine what the likely loss mechanism The electron density in the 
lo torus [Warwick at at., 1979; Bagenal and Sullivan. 1961] decreases 



Ffure 3.1 


The Z>2 intensiUes measured by LET B near closest approach for 
a) Voyager 1 and b) Voyager 2. The Z>2 intensity is piimarily the 
intensity of oxygen and sulfin* ions with energies greater than ~1.1 and 
~0.9 MeV/nuc respectively. Oxygen spectra have been calculated for 
regions indicated in the figure. The number associated with each region 
is the spacecraft -dipole center distance at the time the spacecraft 
crossed the magnetic equatorial plane in that region. 
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sh&rply inside the peak densities at 5.7 inbound (day 64 0924 UT) and 
5.5 Rj outbound (day 64 1420 UT). lo's orbit, on the other hand is at 5.9 
Rj, and the spacecraft crossed the magnetic field lines that pass through 
lo's orbit at 5.65 Rj outbound (day 64 1502 UT). as determined from 
measurements of the lo flux tube [Ness of., 1979a]. 

The energetic oxygen and sulfur losses are therefore more closely 
correlated with plasma features than with lo. Also, it was observed on 
the Voyager 1 inbound pass that the electromagnetic wave activity 
associated with the torus stopped abruptly at ~5.6 Rj (Scarf at., 1979]. 
Hence, the indication is that the losses are due to wave-particle 
interactions in the turbulent region outside 5.6 Rj scattering particles 
into the loss cone. Particle losses and loss mechanisms will be further 
discussed in Chapter 4. Note that mi nima in intensity measurements 
above a fixed energy threshold do not necessarily correspond to minima 
in the phase space densities of the particles. As will be shown (Figure 
3.10), the intensity minima near lo correspond to a change in the radial 
gradient rather than a minimum in the phase space density profile. 

The two decreases at ~9.4 Rj inbound and outbound were caused 
by Europa. In contrast to Voyager 1. the large scale structure in the 
Voyager 2 intensities (Figure 3.1b) is due to varying distance of the 
spacecraft from the magnetic equatorial plane. The peaks at day 190 
1400 UT, 2200-2400 UT, and day 191 0900 UT are at times when the 
spacecraft crossed this plane. 

The inbound Z>2 intensities for the two spacecraft are con^jared 
as a function of radial distance from Jupiter in Figure 3.2. Inside 17 Rj, 
the Voyager 1 and 2 intensities were the same to within a factor of two. 
There is therefore no evidence for major ten^oral variations in the heavy 


Rgure 3J3 


The measured Z>2 intensities lor Voyager 1 and 2 as a function of 
distance from Jupiter. The triangles indicate magnetic equatorial plane 
crossings. 
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ion intensities in this region between the two encounters, and in the 
analysis that follows, both sets of measurements will be considered to be 
of the same relatively stable population. 

3JS Differential Ener^gy Spectra 

Oxygen energy spectra have been calculated from the PHA events 
observed in the ten regions indicated in Figure 3.1 and listed in Table 3.1. 
The 5.3 and 5.9 regions are each the combination of inbound and 
outbound regions symmetrically spaced with respect to the intensity 
minima near lo. For the regions near lo, sodium and sulfur spectra have 
also been calculated. 

Figure 3.3a shows the 5.3 Rj oxygen and sulfur differential energy 
spectra. Also shown is the sodium spectrum from the combined 4.9 and 
5.3 Rj regions. A striking feature in Figure 3.3a is the flattening of the 
oxygen and sulfiir spectra below MeV/nuc and the actual decrease in 
the differential oxygen intensity below 5.5 MeV/nuc. The probable cause 
of this feature is that particle loss mechanisms near lo are more efficient 
for lower energy particles. Further evidence of this effect is provided by 
the observed factor of ~300 decrease in the differential intensity of Zi2 
ions with 1.0 — 1.8 MeV/nuc between ~7 and 5.6 [Kriinigis ef al., 1979]. 
Since the Z>2 integral intensity above ~1.0 MeV/nuc shown in Figure 3.1a 
decreases by only a factor of ~10, it is clear that the low energy ions are 
being preferentially lost in this region. Proton observations between 1 
and 20 MeV near lo [Trainor oL., 1974; Armstrong et al., 1901] also 
indicate preferential low energy losses. 

The 7.9 and 10.1 Rj differential oxygen spectra in Figure 3.3b show 
no evidence of flattening below 7 MeV/nuc. However, trapped energetic 
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Table 3.1 

Region Definitions and Magnetic Field 
Values used in Spectral Analysis 


Region* 

(Rj) 

Voyager 

Time 

(UT) 

Magnetic Equatorial 
Plane Crossing Time 
(UT) 

B Field at 
Crossing 
(Gauss) 

4.9 

1 

64/1106 - 

64/1235 

64/1200 

3.26x10'^ 

5.3 

1 

64/0936 - 

64/1106 

04/1034 

2.79x10“^ 



64/1235 - 

64/1405 

64/1341 


5.9 

1 

64/0842 - 

64/0936 

64/0905 

2.08x10"^ 



64/1405 - 

64/1443 

64/1428 


7.9 

1 

64/0546 - 

64/0659 

64/0608 

7.18x10"^ 

10.1 

2 

190/2214 - 

191/0016 

190/2300 

3.34x10“^ 

12.2 

2 

190/1235 - 

190/1440 

190/1338 

1.74x10“^ 

13.0 

1 

63/2328 - 

64/0110 

64/0008 

1.32x10“^ 

13.7 

2 

190/0907 - 

190/1109 

190/1008 

1.20x10"^ 

16.8 

1 

63/1818 - 

63/2022 

63/1920 

4.85x10"^ 

16.9 

2 

190/0312 - 

190/0517 

190/0415 

6.11x10"^ 

20.0 

2 

- 

f 

189/2238 

3.2 xlO"^ 

20.9 

1 

- 


63/1415 

2.6 xlO"^ 


Distance between the spacecraft and the dipole center at the time of 
the magnetic equatorial plane crossing. Exceptions are 5.3 and 5.9 Rj 
where there were no crossings; tf' listed values are typical distances 
and times in these two regions. 

+ Magnetic field values from N. F. Ness (pnvate communication), 

^ Spectra based on PHA events are not used at 20.0 and 20.9 R, due to the 
difference in Voyager 1 and 2 Z>2 intensities at these radial distances 
and beyond. 
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Figure 3.3 

Differential energy spectra observed in different regions of the 
magnetosphere (Table 3.1). 

a) Oxygen and sulfur spectra from 5.3 Rj and the sodium spectrum from 
the combined 4.9 and 5.3 Rj regions. Note that the sulfur spectrum is 
multiplied by 20. The oxygen and sulfur data were obtained with LET B 
and the sodium data were obtained with LET 6 and D. The lines drawn 
through the oxygen and sulfur data become power-laws of Indices 
approxiriiately -6 at energies greater than ~10 MeV/nuc. The best-fit 
power-law line for the sodium spectrum has an index of -9. The oxygen 
point with dashed error bars at 2.2 MeV/nuc was derived from the 
integral spectrum in Figure 3.6. 

b) Oxygen spectra from 7.9 and 10.1 Rj. The 7.9 data were obtained 
with LET B and U (Voyager 1) and those at 10.1 Rj were obtained with LET 
B (Voyager 2). Power-law lines fit to the spectra have indices of -7. 
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particles do not conserve their enerfiy as they diffuse radially in a 
magnetosphere, so that it is not meaningful to conq)are spectral features 
at a given energy from different radial locations. 

Spectra can be con^ared in terms of other parameters such as 
particle magnetic moments, but before proceeding with such an analysis, 
we will extend the measured differential oxygen spectra to lower energies 
by converting to integral spectra and adding a point at '«! MeV/nuc 
calculated from the Z>2 rate. To do this, it is necessary to determine 
the oxygen fraction of the Z>2 rate. In the 4.9 and 5.3 Rj regions the 
fraction can be directly measured using LI pulse-height information as 
discussed in Chapter 2. The separate oxygen and sulfiu* peaks in Figure 
2.5a can now be understood as the LI energy loss of penetrating 5-6 
MeV/nuc oxygen nuclei (El'-30 MeV) and sulfur nuclei (El~120 MeV). 
This identification is reasonable since the differential energy spectra near 
lo have maxima at these energies. In other regions, the fraction must be 
inferred from elemental abimdance measurements made at higher 
energies, as will now be discussed. 

3.3 blcinenlal Abundances 


At 4.9 and 5.3 Rj, the sulfur to oxygen ratio at constant energy per 
nucleon can be measured with reasonable statistical accuracy between 
approximately 4 and 20 MeV/nuc. The results are shown in Figures 3.4a 
and b. In both regions the ratio is seen to increase toward lower 
energies to values greater than one. The differential energy spectra in 
Figure 3.3a indicate that this strong energy dependence of the ratio is 
caused by particle loss effects near lo. In the other regions listed in 
Table 3.1 there were not enough sulfur events to make a statistically 
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FSCure 3.4 

The energy dependence of the sulfur to oxygen ratio at constant energy 
per nucleon from several r^ions of the magnetosphere. 

a) 4.9 Rj - LET B 

b) 5.3 Rj ' LET B (the same data as shown in Figure 3.3a) 

c) 5.6 - 25. Rj 

Voyager 1 63/0854 - 84/0936 UT LET B and D 
Voyager 2 189/1408 - 191/0410 UT LET B and D 
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si£niAcant determination of the energy dependence of the ratio, so 
instead, all the data from 5.6 to 25 Rj inbound were combined into one 
measurement. The result, shown in panel c, indicates that the sulfur to 
oxygen ratio is approximately independent of energy outside the Z>2 
intensity minima near lo. In order to compare the observed abundance 
ratios to those mode by other instruments at different energies and in 
different locations in the magnetosphere, it is useful to relate them to a 
diffusion theory. 

The radial diffusion theory that will be adopted for the present 
analysis assumes that particles conserve their first and second adiabatic 
invariants but not their third in the diffusion process (see e.g. Schulz and 
Lanzerotti [1974]). Conservation of the first invariant or magnetic 
moment, M = Pj^V2mB where m is the particle mass, B is the magnetic 
field strength and is the momentum perpendicular to the B-field 
direction, gives a linear scaling of energy with magnetic field strength, 
E = MB (nonrelativistic), for the special case of mirroring particles (pitch 
angle. 0=90**). 

Since a particle's magnetic moment is a conserved quantity, 
measurements from different locations in the magnetosphere can be 
compared as a function of M. If analysis is restricted to observations of 
a=90° particles, the magnetic moment can be easily calculated from the 
measured particl>» energies and the local magnetic field strength. 
Abundance measurements such as those shown in Figure 3.4 were not 
restricted in this regard, and, particularly outside 5.0 Rj, represent a 
mixture of pitch angles. However, as will be seen (Figure 3.7), the Z>2 
ions have an equatorial pitch-angle distribution that increases toward 
90®, so that the mixture is biased in favo.r of larger pitch angles. The 
approximation will therefore be made in the analysis that follows that 


46 


M>E/B. The average of the LET B and D pointing directions with respect 
to the magnetic field between 6 and 20 is *^0^ incqplylng that, even ii 
the pitch*angle distribution were istropic in this region, the 
approximation would still be good to ~30%. 

The observed sulfur to oxygen ratio, S/0, is shown as a function of 
M in Figure 3.5a. The error bars with no symbols and with solid symbols 
represent the ratio of the integral sulfur and oxygen measurements 
between and ~15 MeV/nuc. (The exact energy ranges for each point 
are given in the figure caption.). The points with open symbols were 
':.^alculated using published data from the Low Energy Charged Particle 
(LECP) instrument on Voyager 2 (Hamilton ef o£., 1961]. The LECP sulfur 
to oxygen ratios were measured between 10 and 60 Rj inbound for ions 
with energies range 0.60 to 1.15 UeV/nuc. The magnetic moment for 
each measurement (CRS and LECP) was calculated by dividing the lowest 
energy of the measurement range by the local magnetic field value [Ness 
ef of.. 1979a,b and private communication]. Since the sulfur and oxygen 
spectra outside 6 Rj decrease steeply toward higher energies near ~1 

MeV/nuc (-^ocE”® for oxvgen at 10 Rj), the value plotted at a given M is 
dL 

approximately the integral sulfur to oxygen ratio above that M. The 
horizontal error bars on the CRS points represent the variation in the 
magnetic field strength over the measurement region. 

Similar plots of the carbon to oxygen ratio. C/0, and the sodium to 
oxygen ratio, Na/0, are shown in Figures 3.5b and c. The energy 
dependent data inside 5.6 Rj were not used in these figures except for a 
S/0 point for E>6.6 MeV/nuc at 5.3 Rj where the energy dependence is 
less pronoimced, and a Na/0 point for E>6.6 MeV/nuc at 4.9 Rj to make 
use of the few available sodium events. (These are the points at 2.5x10^ 
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Rgure 3.5 

Abundance ratios with respect to oxygen of a) sulfur, b) carbon, and c) 
sodium as functions of the particle magnetic naoments. The data were 
obtained outside 5.6 Rj, except as noted in the text. The magnetic 
moment for each observation was calculated by dividing the lowest 
energy in the measurement energy range by the measured local magnetic 
field strength. The open symbols were calculated iisin£ data from the 
Low Energy Charged Particle (LECP) instrument on Voyager 2 [Hamilton 
ef al.. 1901]. The symbols in the plot are defined as 
Voyager 1 6.6-17.3 MeV/nuc 

• Voyager 2 7.5-11.0 MeV/nuc 

■ Voyager 1 4.7-14.8 MeV/nuc 

♦ = • + ■ 

0 LECP (VGR 2) 0.60-1.15 MeV/nuc. 

The vertical error bars represent statistical uncertainties and the 
horizontal error bars on the CRS data represent the range of magnetic 
field strengths in the region. The power-law lines in the figures are fits 
to the combined CHS and LECP data. The functional forms 8j”e 

|-=1.5M^-^^, ^=2.3xlO-®M°-®®, and ~ = 0 . 1 M ”°-^2 
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MeV/nuc-G in the figure.) 

The agreement between the CRS and LECP obseivations in Figure 

3.5 demonstrates that M is a reasonable parameter to use when 
comparing data from different regions of the the magnetosphere outside 

5.6 The lines shown in the figure are power*law fits to the combined 
data sets. S/0 and Na/0 arc seen to increase towaid smaller magnetic 
moments while C/0 decreases. Therefore an instrument measuring 
abundance ratios in a fixed energy range will detect increasing sulfur and 
sodium and decreasing carbon abundances relative to oxygen as it moves 
toward smaller radial disteuices (larger magnetic fields). Also, an 
instrument at a fixed radial position will measure a sulfur to oxygen and 
sodium to oxygen ratio that increases toward lower energies and a 
carbon to oxygen ratio that decreases, or, in other words, sulfur and 
sodium spectra that are softer than that of oxygen, and a carbon 
spectrum that is harder. If the approximation is made that the energy 
spectra of these elements are proportional to E“^, then the fits shown in 
the figure indicate that 7s = To + 0.2. 7c = 7o“0.6, and 7Ha = ‘yo + 0.1. 

3.4 Integral Energy Spectra 

The data in Figure 3.5 can now be used to calculate the oxygen 
fraction of the Z>2 rate. The calculation '«as performed as follows: The 
LI effective Z>2 thresholds for carbon, oxygen, and sulfur ions in a given 
region were determined from the observed Z>2 rate and pulse-height 
distribution in LI using the information in Table B.2. and were converted 
to magnetic moment thresholds Me, Mq. end Ms, by dividing by the 
observed magnetic field strength in that region (Table 3.1). It was 
assumed that the integral intensity of species k above a given M could be 
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represented by a power law, J]j(>M)=AjjM~^“ where Ajj and 7 ^ are 
constants. The sulfur to oxygen ratio, for exannple. is then 
•^M) = AaM“’^VAoM”^®. With this assurnption, and the observation that, 

inside 20 Rj, carbon, oxygen, and sulfur are more abundant than other 
Z>2 elements, the oxygen fraction of the Z>2 rate is given by 


0+C+S 


M 


JZ>2 rate 


1 + 




IMoJ 


-7o 


|(Ms) 


Mo 


)-7o 


(3.1) 


The carbon to oxygen and sulfur to oxygen ratios were obtained from 
Figure 3.5, and 7 q was determined from the measured spectra. 

The oxygen fraction of the Z>2 rate in each region is listed in 
Table 3.2 in addition to other quantities need to calculate the oxygen 
intensity from the observed Z>2 rate. The oxygen fraction for 4.9 and 
5.3 Rj was not calculated in this way, but was directly measured using LI 
pulse-height information as previously discussed (see e.g. Figure 2.5). 

For each region, an integral oxygen energy spectrum was 
determined, including a point based on the Z>2 rate. Examples from 5.3 
and 10.1 Rj are shown in Figure 3.6. where J(E) is the integral intensity of 
particles with energy greater than E. The higher energy continuous 
portions of the spectra correspond to the differential spectra in Figure 
3.3. and are joined by dashed curves to the points calculated from the 
Z>2 rate. The energy at which the Z>2 point is plotted is the effective LI 
Z>2 threshold for oa^gen ions in the given region, and was determined 
using the data in Table B.2. The value of the effective threshold in each 
region is listed in Table 3.2. 
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Table 3.2 

Quantities Used to Calculate Oxygen 
Intensities from the Z>2 Rate 




Observed 

Effective Oxygen* 

Oxygen Fraction* 

Z>2 Livetime’ 

LZ3 Livetime 

Region 

LET 

Z>2 Rate 

Z>2 Threshold 

of the Z>2 Rate 

Correction 

Correction 

(Rj) 


(S'l) 

(MeV/nuc) 


Factor 

Factor 

4.9 

B 

1.2x10^ 

1.4 

0.52 

1.3 

1.5 

5.3 

B 

1.0x10^ 

1.4 

0.53 

1.2 

1.5 

5.9 

B 

1.3x10^ 

1.5 

0.51 

1.3 

1.5 

7.9 

B 

2.7x10^ 


0.48 

1.7 

1.5 


1.5 




0 

2.5x10^ 



1.6 

1.5 

10.1 

0 

1.9x10^ 

1.3 

0.50 

1.4 

1.6 

12.2 

B 

1.3x10^ 

1.2 

0.50 

1.2 

1.6 


D 

1.1x10^ 



1.2 

1.6 

13.0 

B 

8.0x10^ 

1.1 

0.51 

1.1 

1.6 



1.6 


0 

6.1x10^ 



1.1 

13.7 

B 

7.1x10^ 

1.1 

0.50 

1.1 

1.6 


1.1 

1.6 


D 

6.3x10^ 



16.8 

B 

7.5x10^ 

1.1 

0.53 

1.0 

1.6 





0 

5.4x10*^ 



1.0 

1.6 

16.9 

B 

1.4x10^ 

1.1 

0.52 

1.0 

1.6 




1.6 


0 

1.2x10'’ 



1.0 

20.0 

B 

8.5x10^ 

6.0x10^ 

1.1 

0.53 

1.0 

- 


D 



1.0 

- 

20.9 

B 

1.2x10^ 

7.0x10^ 

1.1 

0.53 

1.0 

1.0 

- 


0 



- 


* Calculated from the observed Z>2 rate and LI pulse-height distribution using Table B.2. 

The nominal oxygen threshold is 1.1 MeV/nuc (sulfur is 0.9 MeV/nuc, carbon is 1.2 MeV/nuc). 
These thresholds include the particle energy loss in the telescope window. 

* From equation (3.1) in all regions except 4.9 and 5.3 Rj, where it was directly measured 
using LI pulse-height information. Note: this is not an abundance fraction. 

’ From Table B.l. 1.0 implies no correction. 

f From Figure B.l. Spectra based on PHA events were not calculated at 20.0 and 20.9 
The 4.9, 5.3, and 5.9 Rj spectra have an additional pulse pileup correction (Appendix 
B) of magnitude 1.4 for Es7 MeV/nuc. 

** The similarity in these numbers is coincidental, due to an increasing carbon fraction 
toward larger radial distances offsetting the effect of a decreasing sulfur fraction. 
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Figure 3.6 

Integral energy spectra of oxygen ions at 5.3 and 10.1 Rj. The 
continuous portions of the spectra are based on PHA events. No binning 
has been done, so that each event is shown as a step up in integral 
intensity. The points with exTor bars were determined from the Z>2 rate. 
The uncertainty in these points is due primarily to abundance 
uncertainties. The dashed line joining the two parts of each spectrum 
represents an interpolation of maximum smoothness [Zmuidzlnas, 1961]. 
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The flattening of the 5.3 Rj spectrum below ~5.5 MeV/nuc 
correspoLds to the decrease in dilTerential intensity seen in Figure 3.3a 
(p. 42). The point with the dashed error bars at 2.2 MeV/nuc in that 
figure was calculated from the integral spectrum in Figure 3.6, and shows 
that the differential intensity decreases even further at energies below 
3.5 MeV/nuc. The 10.1 Rj spect~am continues rising steeply toward lower 
energies between 6.0 and 1.3 MeV/nuc although there is some hardening 
toward lower energies. The spectral index, 7, changes from f at 

6 MeV/nuc to 2 at 1.3 MeV/nuc. The integral spec+ra from all regions will 
be presented and compared in the next section. 


3.5 the Phase Space Density 

3.5.1 Intensity - Phase Space Density Relationships 


Since particle energies change in the diffusion process, the radial 
diffusion theory deals with particle phase space densities at constant 
magnetic moment rather than three-dimensional space densities at 
constant energy. The radial diffusion equation is 


5L - 1 2 ^ 

at ■ ai aL 


+ J(D- X(L) 


(3.2) 


(see e.g. Schulz and Lanzerotti [1974]) where f is the phase space density 
at a particular M, D is the diffusion coefficient, L is the Mcllwain 
parameter, and ^(L) and XiU represent local sources and sinks 
respectively. Equation (3.2) can be used with measurements of the 
gradient of f to obtain quantitative information about diffusive flows and 
particle losses. The objective here will be to determine the radial 
dependence of the energetic oxygen phase space density for several 
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values of M. 

The relationship between the phase space density and the 
measured integral intensity. J. is 


l(p,L.a) 


1 dJ,(E,L) 
d£ 


(3.3) 


where o is the particle pitch angle. Specializing to the a = 90** case and 
changing variables from E and p to M gives 


,(H.L,«=90») = 


(3.4) 


where 7 = - ^ log ^ integral intensity of particles with 

a=90°. The phase space density can therefore be directly calculated 
from measurements of the magnetic field strength, the integral intensity, 
and the logarithmic slope of the integral intensity. 

Equation (3.4) specifies that the integral intensity measurements 
be of mirroring particles. Also, implicit in the diffusion equation is the 
requirement that observations of f at different L values be made at points 
along the dynamic paths that particles follow in the diffusion process. 
This means, for example, that measurements made at radial positions Lj 
of particles mirroring at distances d^ from the magnetic equatoried plane 
should only be compared with one another if the coordinates (Li, di) all 
lie on the mirror point trajectory of the radially diffusing particles. In 
the present study, consideration is restricted to particles mirroring at 
the magnetic equatorial plane to insure that the measurements are along 
miiTor point trajectories, independent of the magnetic field configuration. 

The regions that have been chosen for analysis (Table 3.1) each 
include a crossing of the magnetic equatorial plane, except the 5.3 and 
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I 

5.9 Rj regions, where the spacecraft was near the plane (;S5^ dipole | 

i 

latitude). The exact tinies that the spacecraft crossed the magnetic ; 

equatorial plane were identified using magnetic field magnitude and < 

pointing data (N. F. Ness, private communication) and are listed in Table 

3.1. The L values assigned to each region are the spacecraft -dipole i 

center distance at the time of the crossing, and are the distances used to 

name the regions. 

3.5u3 Pitch Angle Disliibutian Corrections 

The intensity of mirroring particles at a given location can only be 
directly measured with a telescope oriented perpendicular to the local 
magnetic field. At the times of the magnetic equatorial plane crossings, 
however. LET B and D had pointing directions with respect to the field 
that were in general not 90® (see Table 3.3), and, in fact, varied 
considerably from region to region. In order to determine Jj_ from the 
measured intensities, a study was made of the equatorial pitch-angle 
distribution of the energetic oxygen ions. 

It was assumed that the distribution of particles intensities as a 
function of their equatorial pitch angle, oo, could be represented by the 
Legendre polynominal 

f(oo) = ao + a 2 |-(cos2ao - 1-) (3.5) 

where ao and oz constants. This distribution was then integrated with 
the angular response fimctions of the LET B and D Z>2 rate to yield two 
linear equations relating ao end to the measured intensities and 
pointing directions of the telescopes [Gehrels and Chenette, 1981]. The 
ratio a£/ao was calculated from these equations in thirty-minute 



59 


Table 3.3 

Telescope Pointing Directions 
and Correction Factors 


Pointing Angle Z>2 Pointing^ LZ3 Pointing'^' 


Region 

w.r.t. 

Field Line 

Correction Factor 

Correction Factor 

(Rj) 

LET B 

LET D 

LET B 

LET D 

LET B 

LET D 

4.9 

CM 

00 

O 

36° 

1.2 

- n 

1.1 

- f 

5.3 inbound 

60 

35 

1.3 

- 

1.2 

- 

outbound 

65 

32 

1.3 


1.2 

- 

5.9 Inbound 

54 

35 

1.4 

- 

1.4 

- 

outbound 

65 

30 

1.3 

- 

1.2 

- 

7.9 

55 

53 

1.3 

1.4 

1.3 

1.4 

10.1 

65 

46 

1.2 

- 

1.2 

- 

12.2 

77 

41 

1.1 

1.5 

1.0 

1.7 

13.0 

78 

36 

1.1 

1.7 

1.0 

1.9 

13.7 

57 

44 

1.3 

1.5 

1.2 

1.6 

16.8 

49 

50 

1.4 

1.4 

1.4 

1.4 

16.9 

75 

. 53 

1.1 

1.4 

1.0 

1.4 

20.0 

71 

42 

1.1 

1.6 

- 

- 

20.9 

90 

48 

1.1 

1.5 

• 



+ A Compton-Getting correction due to the rotation of the magnetosphere 
Is included in these factors. The approximation was made that the 
rotation speed is 0.8 that of corotation In these regions [McNutt et al., 
1981]. The contribution to the total correction Is 40.05. 

The correction factors for 5.3 and 5.9 Rj take Into account the nonzero 
magnetic latitude of the spacecraft In these regions, as determined from 
the magnetic field data. 

^ In the 4.9, 5.3, 5.9, and 10.1 Rj regions, only LtT B was used because it 
provided adequate statistics by itself and was more nearly perpendicular 
to the magnetic field than LET D. 
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averaginfi intervals from 20 Rj inbound ihroui;h closest approach on 
Voyager 2 and inside 7 Rj on Yoya^er 1. The ratio was found to be 
approximately constant in this region with a distribution of values shown 
in Figure 3.7a. The distribution has a mean value of a^/ao =*-0.7, which 
was chosen as the representative value throughout the region, and a 
standard deviation of 0.2. The shape of a pitch-angle distribution of the 
form given in (3.5) with aa/ao = *~0.7 has a maximum at Oo = 90° as shown 
in Figure 3.7b. Jj. at a given crossing of the magnetic equatorial plane is 
therefore higher than the measured LET B and 0 intensities. 

Using the above pitch-angle distribution and the known telescope 
pointing directions (Table 3.3), the value of at each magnetic 
equatorial plane crossing was calculated from the LET 5 and D intensities. 
The ratio between and the measured intensities is called the pointing 
correction factor, and is listed in Table 3.3. The pointing correction 
factors for the Z>2 rate are slightly different than those for the LZ3 rate 
due to their different angular response functions. 

The specific procedure followed in dete rmining Jj,(M) in each region 
was as follows; Integral energv spectra were calculated for LET B and D 
(in some regions only LET B) using the values of the Z>2 rate and LZ3 
rate (equation 2.1) measured at the magnetic equatorial crossing. The 
pointing corrections (Table 3.3) were then applied, and the LET B and D 
spectra combined, to 3 rield a single spectnun of 00=90** oxygen ions in 
each region. In the 4.9, 5.3, 5.9, and 10.1 Rj regions, only LET B was used 
to determine because it provided adequate statistics by itself in these 
regions and was pointed more nearly perpendicular to the magnetic field 
direction than LET D (Table 3.3). The value of the magnetic field strength 
used to convert the integral energy spectrum in each region into an 
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Rgure 3.7 

ft) A histogram of measurements of the pitch^ani^le distribution 
parameter ratio a^^/aQ. The distribution of particles intensities as a 
function of their equatorial pitch angle, oq. is assumed to be of the form 

f(oo) = ao + aa^(cos*ao-^) The data shown in the figure represent 

thirty-minute averages of a^/ao as determined by the LET B and D Z>2 
rates measured by Voyager 2 on day 190 and by Voyager I between 
64/0730 XJT and 64/1630 UT. Periods in which the telescope pointing 
directions relative to the magnetic field direction were within 10® of one 
another were not included due to the large uncertainties in the 
calculated ratios during them. The mean value of the ratio for the data 
shown in the figure is -0.67. It was therefore assumed that a2/ao=~0.7 
throughout this region. 

b) The shape of the pitch-angle distribution for a2/€^)=-0.7. The ratio, 
ae/fto, was directly measured only for 30®i ao^lSO®, so that the shape of 
the actual distribution near Oq = 0® and 180® is not necessarily as shown. 



Number of Intervals 
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integral magnetic moment spectrum was the measured value at the time 
of the magnetic equatorial plane crossing and is listed in Table 3.1. 


3.^.3 Phase ^ace Density Deternanation 


The phase space density of the energetic oxygen ions can now be 
determined from the measurements of Jj^(M.L) and B(L) using equation 
(3.4) 




As a first step, we show in Figure 3.8 the quantity as a function of 

magnetic moment for the measured spectra in all regions. The event 
portions of the Voyager 1 13.0 and 16.0 Rj spectra are represented by 
single points instead of by complete spectra due to the few PHA events 
obtained in these regions. These measurements, as well as the two 
intensity measurements based on the Z>2 rate, are ingiortant, however, 
becaixse they confirm that, to within a factor of two, the Voyager 1 and 2 
intensities of energetic oxygen ions were the same between 12 and 17 Rj. 
Therefore, inside 17 Rj. spectra from the two spacecraft will be directly 
compared in the analysis that follows. Outside 17 Rj, the Voyager 1 
intensities of energetic Z>2 ions were higher than those of Voyager 2, as 
can be seen by the two intensity measurements neai’ 20 Rj in this figure 
and by the radial con^jarison of the Z>2 intensities in Figure 3.2. 

The short-dashed curves in Figure 3.0 are extrapolations of the 
measured spectra to smaller magnetic nx)n,r!nts, while the error bars 
represent the range of uncertainty in the extrapolations. Except at 10.1 
Rj, the upper and lower limits for each error bsu- were determined by 


Figure 3.B 


Integral oxygen spectra from all regions listed in Table 3.1 divided by the 
square of the local magnetic field strengths, as a function of magnetic 
moment. The spectra are represented in the same way as those in Figure 
3.6. The event portion of the 4.9 Rj spectrum is shown as a dotted curve 
to distinguish it from the 5.3 Rj spectrum Intensities derived from the 
PHA events at 13.0 and 16.8 Rj are represented by open symbols rather 
than by line spectra because few events were obtained in these regions 
(''i'lO events in each region), and are shown primarily for comparison with 
the Voyager 2 spectra. The short-dasxhed curves and the dash-dotted 
curve at 10.1 Rj represent extrapolations of the measured spectra as 
discussed in the text. The error bars without s 3 Tnbols indicate the 
uncertainties in the extrapolations. 
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extending the spectra from the Z>2 point with lines of slope 72>s and 

1-^ respectively, where 7z>z is the slope at the Z>2 point. 

At 10.1 Rj. a longer extrapolation was performed than in other 
regions in order to extend the spectrum to 70 MeV/nuc-G (see section 
4.6). In this case, the shape of the spectrum was calc:ulated using the 
energy per charge spectrum measured at these magnetic nooments by 
LECP at ~16 Rj on Voyager 2 [Hamilton at al., 1961]. Hamilton at at. 
[1861] show that different elements have similar spectral shapes at the 
same energy per charge, if one assumes reasonable charge states for tiie 
different ions. They therefore combine all elements into a single 
spectrum that covers approximately two orders of magnitude in energy. 
The lowest energy oxygen measurement in the spectrum is at M«500 
MeV/nuc-G, but the spectral shape is determined by protons and helium 
ions to magnetic moments as low as ~60 MeV/nuc-G. The integral oxygen 
intensity profile used to extrapolate the 10.1 Rj spectrum was obtained 
by converting from energy per charge to energy per nucleon, integrating 
to obtain J(E), and then converting to J(M) using Ba^exlO"^ G. 

LECP observations were also used as a consistency check of the 
spectral slopes at the Z>2 points in Figure 3.6. Figure 3.9 shows these 
slopes compared to those measured in the energy range from 0.60 to 3.17 
MeV/nuc between 10 and 25 Rj [Hamilton at ol., 196 i]. The values from 
the two instruments are the same within statistics, and both show the 
same trend of softer spectra toward increasing radial distance. 

In order to use the diffusion equation (3.2) to understand particle 
observations, one must determine the radial dependence of the phase 
space density at a constant value of M. Radial profiles of f at constant 
magnetic moment can be directly obtained for the energetic oxygen ions 
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Rgure 3.9 

A comparison of the slope of the oxygen spectrum at the Z>2 point (~1 
MeV/nuc) observed by CRS, with the power-law index between 0.60 and 
3.17 MeV/nuc measured by LECP [Hamilton 9t al., 19811. Note that the 7 

dJ 

plotted is the power-law index for a differential spectrum so 

that, for the CRS observations, 7 = 72>a+l. where 72>2 the slope of the 
Integral spectrum at the Z>2 point. 
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from Figure 3.8 using equation (3.4). The results are shown for several 
values of U in Figure 3.10. The filled S 3 rmbols in the figure denote phase 
space densities detennined directly from the spectra in Figure 3.6, while 
the open symbols are based on the extrapolations of the spectra 
discussed above. 

The uncertainties, shown for each measurement of f, represent 
xmcertainties in the magnetic field strength, in the oxygen 

abundance fraction of the Z>2 rate, a,i,und«nce* pulse pileup 

correction factor for the two-detector portions of the spectra near lo, 

in 7(M,L), 07 . in the pointing correction factor, and 

statistical counting uncertainties in the number of events in the event 
portions of the spectra, aerenis’ The estimated contribution to Of from 
each of these sources is listed in Table 3.4 for all measurements. An 
additionaJ '«15% systematic uncertainty in the livetime correction factor 
applies to all measurements. 

The observations at all values of M show a positive radial gradient 
in the phase space density of the energetic oxygen ions between 5.6 and 
17 Rj. The diffusive flow in this region is therefore inward. Inside 5.8 Rj, 
the phase space density is nearly constant. The discontinuity in the 
gradient of f at ~5.6 Rj corresponds to the minimn in the Z>2 intensities 
near lo (Figure 3.1). 
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Figure 3.10 

The radial dependence of the phase space density of the energetic 
oxygen ions at constant magnetic moment, for several values of M. The 
filled points are derived from the measured portions of the spectra shown 
in Figure 3.8, while the open points are based on the extrapolations. The 
4.9 and 5.3 Rj measurements at 6.9x10^ MeV/nuc-G were averaged 
together to improve statistics. Solid lines join the measurements at given 
magnetic moments between 7.9 and 16.9 Rj. The 5.3, 5.9, and 7.9 Rj 
measurements are connected by dashed lines to indicate that the radial 
profile in this region may have Oiscontinuities between the points. 
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Table 3.4 

Contributions to the Phase 
Space Density Uncertainties 


M 

(MeV/nuc-G) 

L 

"l/B* 

(t) 

°abund. 

(«) 

"events 

(t) 

"2-0 

(») 

a 

r 

M 

"point. 

{*) 

f 

(relative) 

"f 

(!) 

7.0x10^ 

4.9 

6 

15 

0 

0 

15 

10 

1.4x10'^ 

24 


5.3 

« 

N 

« 

M 

U 

M 

1.7x10’^ 

M 


5.9 

M 

H 

M 

N 

M 

M 

1.5x10° 

H 

2.0x10^ 

4.9 

6 

0 

7 

15 

15 

10 

7.5x10"^ 

25 


5.3 

M 

M 

6 


N 

M 

9.8x10"^ 

M 


5.9 

n 

m 

7 

M 

■ 

M 

1.3x10“^ 

M 


7.9 

« 

15 

0 

0 

M 

« 

9.1x10° 

24 


10.1 

- 

- 

extrapolation 

- 

- 

1.3x10^ 

36 


12.2 

- 

- 

extrapolation 

- 

- 

6.3x10^ 

72 

4.0x10^ 

4.9 

6 

0 

27 

0 

20 

10 

3.3x10’^ 

36 


5.3 


M 

21 

M 

H 

M 

3.9x10’^ 

31 


5.9 

M 

M 

30 

II 

« 

II 

9.4x10"^ 

38 


7.9 

II 

15 

0 

M 

15 

M 

8.4x10"^ 

24 


10.1 

m 

H 

M 

M 

H 

II 

7.7x10° , 

U 


12.2 

- 

- 

extrapolation 

- 

- 

4.1x10^ 

40 


13.7 

- 

- 

extrapolation 

- 

- 

8.0x10^ 

54 

6,9x10^ <4.9,5.3> 

6 

0 

100 

0 

20 

10 

7.8x10"^ 

103 


5.9 

n 

M 

M 

N 

a 

« 

3.6x10'^ 

M 


7.9 

« 

■ 15 

0 

N 

15 

M 

8.7x10'^ 

24 


10.1 

M 

M 

N 

M 

M 

M 

8.8x10"^ 

M 


12.2 

M 

M 

M 

N 

H 

M 

5.7x10° 

M 


13.7 

- 

• 

extrapolation 

- 

- 

9.6x10° 

30 


16.9 

- 

- 

extrapolation 

- 

- 

3.5x10^ 

67 

9.4x10^ 

5.9 

6 

0 

150 

0 

50 

10 

1.3x10'^ 

160 


7.9 

H 

k 

18 

II 

20 

II 

1.6x10'^ 

29 


10.1 

• 

15 

0 

0 

15 

■ 

2.6x10*^ 

24 


12.2 

• 

m 

II 

• 

II 


1.7x10° 

m 


13.7 

M 


» 

II 

II 

II 

3.5x10° 

H 


16.9 

- 

- 

extrapolation 

- 

- 

9.3x10° 

49 

1.7x10^ 

7.9 

6 

0 

300 

0 

50 

10 

4.7x10"^ 

304 


10.1 

■ 

a 

21 

II 

15 

II 

1.3x10’^ 

28 


12.2 

■ 

15 

0 

II 

• 

• 

1.4x10*^ 

24 


13.7 

■ 

• 

II 

«■ 

H 

n 

3.1x10'^ 

« 


16.9 

M 

M 

II 

II 

•1 

H 

8.1x10"^ 

M 
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Chapter 4 
Discussion 


4.1 Ove r vi ew 

The measurements of the phase space density of the energetic 
oxygen ions (Figure 3.10) can be used both to learn more about the ions 
themselves and to place limits on parameters that control the energy and 
particle budget of the magnetosphere. Concerning the ions, the specific 
objectives of this chapter are to 1) establish their overall fioir pattern 
and the location of their acceleration region. 2) calculate the total 
inward diffusive flow between 5 and 1 7 Rj. and 3) determine the 
magnitude and cause of the losses in the region. Objectives 2 and 3 
require knowledge of the diffusion coefficient. The energetic oxygen data 
can bfe used to place limits or this parameter that are adequate for the 
analysis, and that are also interesting in that they specify limits on the 
total plasma outfiow from the torus and thereby on the rate at which 
oxygen and sulfur ions are produced near lo. 

The radial dependence of the inward diffusive flow rate of the 
oxygen ions, calculated using the diffusion coefficient limits, indicates 
that almost all of the ions are lost between 12 and 5 Rj. In order to 
determine the cause of these losses, we calculate in the fifth section of 
the chapter the relative in^ortance of several possible loss processes. 
The only viable loss mechanism is found to be pitch-angle scattering into 
the loss cone. Tliis iir^lies that most of the energy delivered to the inner 
magnetosphere by the oxygen and sulfur ions ends up in the Jovian 
atmosphere. It is of particular interest to compare this power input into 
the atmosphere with the total power required to produce the observed 
auroral activity [Broadfoot «f ai., 1979; Sandel at ai., 1979], which is the 
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most enersetic phenomenon associated viih the inner magnetosphere. 

4J2 Hie 0:^gen Flow Pattern and Acceleration Prooen 

The positive radial gradient in the phase space density of the 
energetic oxygen ions between 6 and 17 Rj causes an inward diffusive flow 
in this region, and indicates that the acceleration region is outside 17 Rj. 
Oxygen ions with magnetic moments in the range 1C*~10* McV/nuc-G 
have also been observed throughout the magnetosphere by the I£CP 
instrument [Krimigis al., 1979; Hamilton et al,, 1961; Krimigis *t al., 
1981]. In addition to the *^16 Rj measurement discussed in the last 
chapter, LECP observations in two other regions are of particular interest 
in this study; one was made just inside the magnetopause at 65 Rj on 
Voyager 1 [Krimigis et al., 1979] and the other in the magnetospheric 
wind region (r>150 Rj outbotmd) on both Voyager 1 and 2 [Krimigis et al., 
1981]. The J/B^ calculated from the 65 Rj measurement is two to three 
orders of magnitude larger, at the same value of M, than that at 16.9 Rr 
in Figure 3.8. This observation does not necessarily imply particle flow 
directly in from 65 to 17 Rj since there may be a maximum in the phase 
space density at some intermediate distance, but it does confirm that the 
acceleration region is outside 17 Rj. The magnetospheric wind 
observation shows that ions flow outward as well as inward from the 
acceleration region. 

The oxygen and sulfur rich composition of the energetic ions 
indicates that their original source is the lo torus. The observations 
therefore suggest that plasma ions diffuse outward from the lo torus, are 
nonadiabatically accelerated in some region outside 17 R;, and then 
diffuse inward and outward from the acceleration region. The radial 
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diffusion process will change particle ^luergics. but is assumed to 
conserve magnetic moments. The LECP measurements just described and 
Uie CRS measurements reported here both require an acceleration 
mechanism that is capable of increasing particle magnetic moments from 
that of the plasma near lo, McV/nuc-G, up to *^10* MeV/nuc-G. 

One possible mechanism is that proposed by Eviatar et al. [ 1976] to 
predict the observed energetic sodium component in the magnetosphere. 
In that model, fast neutral particles, possibly produced by charge 
exchange between neutrals and corotating plasma near lo. escape the 
magnetosphere and become ionized in the solar wind. The ions are then 
accelerated to the ~400 km/ s soleur wind speed and re-enter the 
magnetosphere. This mechanism will, however, result in magnetic 
moments of not more tiian ~20 MeV/nuc— G. Two other possibilities arc 
stochastic acceleration by cyclotron waves in the plasma of the outer 
magnetosphere [Papadopoulos et ai., 1960], or the application ol 
magnetic punq)ing models [Sentman et al., 1975; Goertz, 1976; Borovsky 
sf 04., 1981] to heavy ions. 

4.3 Diffusion Coefficient Deternsnations 

In order to use the observed oxygen spectra and phase space 
density profiles to determine the diffusive flow rate at which particles and 
energy enter the inner magnetosphere, it is necessary to know the 
diffusion coefficient The objective of this section is to place limits on the 
value of the diffusion ccefficient between 5.3 and 17 Rj based on limitmg 
assumptions concerning the magnitude of the losses in the region. 

The diffusion equation (3.2) for a steady-state condition in a region 
in which there are no sources and where the loss term is of the form 
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T, where r ia ihe particle lifetime, is 



L« aL 



(4.1) 


We first consider the extreme case in which there are no losses (r-*<») 
between 5 and 17 Rj. 


4.S.1 Lossless DilfusiDn 


For the case of no losses, equation (4.1) gives 


D 


cL^ 

df/aL 


(4.2) 


where c is a constant. Relative diffusion coefficients can therefore be 
directly calculated from the phase space density measurements. For 
ease of calculation, the radial dependence of f at a given magnetic 
moment is approximated as a series of power-law lines connecting the 
measured points at different values of L. Coirporison with the smoothed 
curves in Figure 3.10 indicates that this is a reasonabl'.’: approximation. 

Figure 4.1 shows the values of the diffusion coefficient at five radial 
positions calculated using equation (4.2). The points are normalized to 
D=1.0 at L=11.4, and have errors that represent the uncertainties in the 
measured phase space densities. The relative diffusion coefficients for 
different magnetic moments are seen to be the same within statistics. 
Also shown in the figure arc two lines of the form D=DoL°, where Dq and n 
are constants. The radial dependence of the calculated points can be 
adequately represented by n=— 7 between *^7 and *^12 Rj and n=0 between 
~12 and *^15 Rj. This analysis shows that it is possible to obtain radial 
profiles of f between 7 and 15 Rj similar to those measured, without 
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Figure 4.1 

Relative values of the lossless diffusion coefficient between 6 and 17 Rj at 
several oxygen magnetic moments: calculated using equation (4J3). The 
gradient in the phase space density was determined from the measured 
values of f (Figure 3.10) assuming f«L*t between measurements. The 
points are normalized to D=1.0 at L=11.4, and have errors that represent 
the uncertainties in f. Power-law lines of the form D = DqL“ with n=— 7 
and n=0 are shown for conparison. 
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invoking any particle losses- The resulting radial dependence of D, 
however, does not satisfy a desired constraint that will now be discussed 
The radial dependence of the diffusion coefficient in the inner 
nsagnetosphere can be independently determined using Voyager 
measurements of the plasma density profile [Froidevaux, 1960; 
Richardson at al., 1930. Richardson and Siscoe. 1981; Siscoe at at, 1981]. 
In particular, the recent analysis by Siscoe at al. [1961] between 5.8 and 
9 Rj gives 

D = li.8xl0-®® S 5.8<Lsi7.12 

D = 1.7x10-^® S s~^ 7.13<L^8 (4.3) 

D = 2.6x10-®’ S s-^ 8<Ls9 

where S (in ions/s) is the lo source strength of oxygen and sulfur ions. 
As discussed in the introduction, S is probably in the range from 10^ to 
10“ ions/s, but may be as large as 10^ ions/s. The most likely diffusion 
mechanism between 6 and 9 Rj is flux tube interchange driven by 
centrifugal forces on the plasma created near lo [Siscce and Summei's, 
1981], and therefore, since the mechanism involves magnetic field motion, 
the diffusion coefficient is the same for the plasma and the energetic 
ions. 

The diffusion coefficient in (4.3) is not adequate by itself for the 
present analysis bece.ise of its limited radial range (L<9), and because of 
the uncertainty in S, but it does indicate that the n=-7 radial 
dependence of the lossless diffusion coefficient between ~7 and ~12 Rj in 
Figure 4.1 is inconsistent with the plasma measurements. Therefore, 
diffusion with non'^zero losses will now be considered. 
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4.3.2 Lossy niffusioii 

As &n opposite limit to the cese of no losses, it is assumed in this 
section that strong pitch-angle diffusion is occurring. This means that 
there is sufficient pitch-angle scattering to refill the loss cone with 
particles as quickly as it empties, and therefore that the lifetime of the 
particles to pitch-angle scattering is a minimum. The loss cone solid 
angle at the magnetic equatorial plane (both directions) is approximately 
Swaoj* where aoj is the equatorial loss cone angle, and the time required 
to empty it of particles is approximately one quarter of tlie bounce 
period, Tb, [Kennel emd Petschek, 1966], so that 

7V^ 4 

T«t«>ng = • (■^•4) 

OtQl ' 2 

For nonrelativistic particles at constant M, assuming a dipolar Jovian 
magnetic field of 4 G-Rj®, 

« 5 ir'''* [4-3/L]*'* s (4,6) 

where [M]=MeV/nuc-G. The deviation of the measured magnetic field 
from a dipole [Connemey et ai., 1901] tends to decrease oq] and thereby 
increase The difference, however, is estimated to be less tlian 

~50% between 6 and 17 Rj, and in the analysis that follows, equation (4.5) 
will be used to calculate 

The objective now is to calculate the diffusion coefficient using 
Inside 9 Rj, the radial dependence of D is specified by Voyager 
plasma measrirements and is given explicitly by equation (4.3). In this 
region therefore, the diffusion coefficient can be fully determined by 
obtaining a value for S. Between 7 and 0 Rj there are two discontinuities 
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in the plasma diffusion coefficient representinj; observed discontinu'i js 
the radial gradient of the plasma density [Siscoe et al., 1961]. T'.ih 
implies that there eire also discontinuities in the radial gradient of the 
energetic oxygen phase space density between the 5.9 and 7.9 R; 
measurements, and therefore that the radial resolution of the energetic 
oxygen measurements is insufficient in this region to independently 
determine the gradient in f for a comparison with the plasma 
obsci^'^’ons. Between 8 and 9 Rj, however, a comparison can be made 
and a value obtained for the diffusion coefficient. Outside 9 Rj, a 
numerical solution of the diffusion equation relates the diffusion 
coefficient to the one calculated between 6 and 9 Rj. 

The radial dependence of f in this analysis is again approxinwted 
as a series of power-law lines connecting the measured points in Figure 
3.10. Specifically, assuming f«L*i and D=2.6xlO“®^SL^ between 8 and 9 Rj, 
the diffusion coefficient at 8.9 is 


D(8.9) 


8jf 


(4.6) 


For a given M, the values of the diffusion coefficient at larger radial 
distances were determined from the value at 8.9 Rj using a finite 
difference solution of the diffusion equation 14. l). For instance, D(11.4 
Rj) was calculated using 


D(1 1.4) |f-{l 1.4)/ 11.42 _ :)(8.9)|f^8.9)/8.92 


3f 


lo.r 


aL 


dL 


11.4 “ 8.9 


] 


^ tdO.Q 


(4.7) 


Figure 4.2 shews the values of D obtained in this way at several 



Figure 4J3 


Values of the lossy dilTusion coefficient between 9 and 17 Rj for the case 
of strong pitch-angle diffusion (T=T,troiia: given by equation (4.5)), 

The points at 8.9 Rj were calculated using equation (4.6), where q is the 
power-law index of the phase space density between 7.9 and 10.1 Rj. This 
equation was derived assiuning that D«L* between 6 and 9 Rj as 
indicated by Voyager plasma data [Siscoe et al., 1981]. The solid lines in 
the figure show the dependence of the plasma diffusion coefficient on the 
lo source strength, S, of oxygen and sulfur ions. The points at 11.4. 13.0, 
and 14.4 were calculated using equations like (4.7). The dashed lines are 
the diffusion coefficient limits (Dy - upper, D] - lower) that will be used in 
the analysis that follows, with the range of values at 5.3 Rj determined 
from Pioneer measurements of energetic protons and electrons [Goertz cf 
al., 1978; Thomsen et al., 1977]. 
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radial positions for different values of M. The diffusion coefficient at a 
given radial position and the general radial dependence of D are seen to 
be the same within statistics for all magnetic moments. It is important to 
note that, since the diffusion coefficient was calculated using Tgtroqgi it 
represents an upper limit (assuming losses are caused by pitch^angle 
scattering). Also shown in the figure is the plasma diffusion coefficient of 
equation (4.3) between 8 and 9 Rj for S= 10*’', 10*®, and 10** icns/s. A 
comparison of these lines and the upper limits obtained above indicate 
that the lo source strength of oxygen and sulfur ions is less than 
•^10*® ions/s. 

As discussed in Chapter 1. S is most likely greater than 10*^ ions/s, 
which places a lower limit on the diffusion coefficient of approximately 
one tenth the value of the upper limit. This inyjlies that the lifetime of 
the energetic oxygen ions between 8 and 17 Rj is within em order of 
magnitude of r,tron«. The range of diffusion coefficient that will be used 
in further analysis is represented by the dashed lines in Figure 4.2. The 
limits shown at 5.3 Rj were not determined from the energetic oxygen 
data, but represent the range of values determined from energetic proton 
and electron measurements (Goertz et oL., 1979; Thomsen et al., 1977). 

4.4 Diffusive Flow Rates 

The diffusion coefficient limits shown in Figure 4.2 can now be used 
with measurements of the energetic oxygen number density and phase 
space density gradient to determine the inward flow rate. The diffusive 
flow, I, of ions with magnetic moment greater than M across radius r is 
given by 
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I>M = n>||-27rr'2h-voRj^ ions/ s (4.8) 

where n>ii is the number density (cm‘*^ of oxygen ions with magnetic 
moment greater than M calculated from the measured spectra, h is the 
scale height (Hj) of the particle distribution relative to the magnetic 
equatorial plane, and Vq is the average radial diffusion velocity (Rj/s) 

given by vd = D*(-^-|~) (see e.g. Schulz and Lanzerotti [1974]). 

In a dipolar magnetic field, the value of h for a particle population 
with a constant pitch-angle distribution depends on radial distance, 
increasing with increasing r. The scale height calculated for the 
energetic oxygen pitch-angle distribution (Figure 3.7b) using the actual 
magnetic field measurements was found, however, to be approximately 
constant between 6 and 17 Rj, equal to '-3 Rj. This is due to the 
distortion of the magnetic field at larger radial distances caused by 
currents that flow in the plasma disk near the magnetic equatorial plane. 
The magnetic field strength increases more rapidly as one moves away 
from the plane than for a dipole field, ana, therefore, at larger r, the 
scale height of the energetic particles is reduced. 

Figure 4.3b shows values of I>n obtained from equation (4.8) for 
oxygen ions with magnetic moment greater than 4.0x10^ and 1.7x10® 
MeV/nuc-G between 5.3 and 17 Rj. The error bars represent all 

1 df 

uncertainties in Table 3.4 (except oy), a 15% uncertainty in ^ 

50% uncertainty in the relative values of D at different radial locations. 

1 A# 

(Larger uncertainties in D and were used at 5.3 Rj.) An additional 

factor of 3 absolute uncertainty in D applies to ail measurements. The 
charge density profile of the plasma in this region, measured by the 
plasma instrument on Voyager 1 [Bagenal and Sullivan, 1981], is shown 
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Figure 4.3 

a) The charge density of the plasma in the inner magnetosphere as 
measured by the Voyager plasma instrument [Bagenal and Sullivan, 1981]. 

b) The inward flow rate of energeUc oxygen ions of magnetic moment 
greater than 4.0x10* and 1.7x10* MeV/nuc-G calculated using equation 
(4.6). In addition to the uncertainties shown in the flgure (discussed in 
the text) there is an absolute uncertainty of a factor of ~3 that applies 
to all measurements. 
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for consparison in panel a and will be discussed in section 4.5. L 

It is evident from the fl£ure that the total number of energetic 
oxygen ions flowing into the inner magnetosphere is small ions/s 

at 10 Rj for M=4.0xl0® MeV/nuc-G) compared with lo’s source strength of 
ions/s. Another observation is that almost all of the ions 
flowing into this region are lost between 6 and 12 Rj. The region with the 
largest losses in tei^ms of total number of particles is 8 to 12 Rj. The 
profile of the inward flow rate for M>4.0xl0* MeV/nuc-G ions between 5.3 
euid 5.6 Rj is drawn approximately flat to indicate that the losses are 
most likely much less in this region than between 5.6 and 7.9 Rj. 

An estimate of the average lifetime, t, between 5.6 and 7,9 Rj can 
be obtained from the measured values of I using 

T= -t/ log* |[(5.6)/I(7.9)|, where t is the time spent diffusing from 5.6 to 

7.9 Rj. For D«10“® s~* and •^■|7-«2, one finds t«1.2xl0® s and 

f oL 

t«1.4x10® s, to be compared with t^,.o.t.,(6.5Rj)= 1.4x10* s. The average 
lifetime is therefore approximately an order of magnitude longer than the 
strong pitch'angle scattering lifetime in this region. Note that, even 
though there are relatively few particles lost outside 12 Rj, the lifetime is 
still close to in this region. For instance, repeating the above 

calculation between 12.2 and 13.7 Rj for M>1.7xl0® MeV/nuc-G ions gives 
7^1.7x10® s compared with Trtrom(13Rj)=3.5xlO® s. The reason fewer 
particles are lost in this region is that the ions are diffusing more rapidly 
(D is larger), the loss cone is smaller, and the bounce time is longer than 
at smaller radial distances. 

Although the energetic ions that flow into the inner magnetosphere 
are small in number relative to the number of ions produced near lo, 
they do bring a significant amoimt of energy into the region as will be 
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discussed in section 4.6. It is therefore important to understand what 
the loss mechanism is in order to determine where the energy is 
delivered. 

4.5 Particle Lon Mechanism 

It has been implicitly assumed in several previous sections of this 
thesis that the predominant loss mechanism for the energetic oxygen 
ions inside 20 Rj is pitch^angle scattering into the loss cone. In this 
section, the lifetime for several other possible loss mechanisms is 
estimated to see if they are significant in this region. The particular 
processes considered are geometric absorption by lo. absorption by 
micron-sized dust grains, and energy loss in the plasma of the inner 
magnetosphere. Energy-loss processes produce losses by shifting particle 
energies below analysis thresholds. They also represent real losses in the 
sense that energy is transferred from the particle to the medium In 
order to calculate the energy loss in the plasma, information concerning 
the charge state of the energetic ions is needed, and therefore, charge 
state limits will also be determined. The observed losses indicate that the 
effective lifetime for the energetic ions is at most ten times the strong 
pitch-angle diffusion lifetime, so that calculated lifetimes for 1 MeV/nuc 
ions should be compared with ~9 days at 0 Rj (10 t«,^ for 128 
MeV/nuc-G particles) and ~3 days at 6 Rj for 54 MeV/nuc-G 

particles). 

4.5.1 Geometric Absorption by lo 

In section 3.1 it was seen that the mmima in the Z>2 rate near lo 
occxirred at 5.6 Rj inbound and 5.4 Rj outbound, well inside lo’s orbital 
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position ol 5.9 The locations of the minima, which specify the inner 
edge of the lossy region, were found to be correlated with the locations 
of discontinuities in plasma density and wave activity. Figure 4.3 gives 
further evidence that the losses are not caused by lo. The total number 
of inflowing ions begins to decrease already at 12 Rj and falls steeply 
between 10 and 6 Rj. The plasma charge density proflle, shown in panel 
a. indicates that the region of high plasma density extends well beyond 
the orbit of lo, to Rj. Also, significant wave activity was observed by 
the plasma wave instruinent on Voyager 1 and 2 [Scarf et of., 1979; 
Gurnett at ai., 1979; Scarf at at., 1961] out to at least 20 Rj, with 
particularly Lxtense activity between 5.6 and 1C Rj. It is therefore found 
once again that the lossy region is better correlated with regions ol 
enhanced wove activity and plasma density than with the orbit of lo, and 
that certainly, since most of the energetic ions are lost between 12 and 6 
Rj, Io does not absorb a significant fraction of the total oxygen energy 
flow into the inner magnetosphere. 

4.fl>2 Absorption by Dust 

It is postulated that SOj dust grains of '-O.l /un radius escape from 
lo's volcanic plumes into the Jovian magnetosphere [Johnson at ai., I960]. 
The number density of the grains near lo is estimated to be 
n^j^«10“® cm*^ [Morfill at ai., I960]. Therefore, the collision JileUme, t. 
of a particle of velocity v (=1.4x10® cm/s lor E=1 MeV/nuc) m this 
regions is 

T = at 10 years (4.9) 


where a is the cross-sectional area o! a dust grain (~3xl0"*° cn:?). A 1 
MeV/nuc oxygen ion can penetrate approximately ten grains before losing 
too keV/nuc, but the life Jme fox' collision with even one grain is so long 
that it is clear absorption by dust is an insignificeuit loss mechanism. 

4.5.3 Oxygen Charge State 

An estimate of the charge state of the energetic oxygen ions is 
needed to calculate the rate at which the ions lose energy in a plasma 
(section 4.5.4). the esUmate is made assuming that the only process 
altering the charge state of the oxygen ions as they diffuse through the 
inner magnetosphere is electron-impact ionization due to free plasma 
electrons. Since ionization due to plasma ions and neutral gas in the 
region is not included, the calculated charge state represents a bwer 
limit. It is also a lower limit because the assunq>tbn is made that the 
ions are singly ionized when they cross 10 Rj. If the energetic ions 
outside 10 Rj have the s ame charge state composition as the plasma ions 
in the hot torus, then they are a mix of O’** and [Bagenal and 
Sulhvan, 1981]. 

An approximation made in the calculation is that the electron- 
capture cross sections are negligible compared with the ionization cross 
sections. The average charge state that an ion with velocity fi (=v/c) 
acquires after imdergoing a large number of interactions with a medium 
is given by 

Z* = z[l -exp( -125/r/ Z^^)\ (4.10) 

[Barkas, 1963] where Z is the nuclear charge, and represents the charge 
state at which ionization and electron-capture cross sections are equal. 
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Therefore, since Z*>6 for >1 MeV/nuc oxygen ions, electron capture can 
certainly be neglected for ions of charge state less than ~5. 

Hie ionization lifetimes were calculated as a function of radius 
between 5.6 and 10 Rj using experimental and theoretical electron-impact 
ionization cross sections [Lotz. 1967] and electron densities inferred from 
plasma wave observations [Gumett cf at.. 1961]. The measured diffusion 
coefficient and phase space density gradient were used to determine the 
length of time the energetic ions spend at each radial distance. The 
resulting ionization state lower limit for 400 MeV/nuc-G oxygen ions is 
shown in Figure 4.4a as a function of position. Since the calculation was 
performed using average ionization lifetimes, the figure represents the 
average position at which a given ionization state occurs. The results 
indicate that the plasma is dense enough in the inner magnetosphere to 
strip the energetic ions of most of their electrons. 

The oxygen gyrofrequency. which is a parameter of importance fcr 
studies of wave-particle interactions, is shown in Figure 4.4b for 400 
MeV/nuc-G ions with charge states given by panel a. Particle 
gyrofrequencies are proportional to QB. where Q is their charge state, 
and therefore the increase seen in the figure between 10 and 5 Rj is due 
both to increasing magnetic field strength and charge state. 

4b5.4 Snefgy Lom in Plasm 

Figure 4.4a shows that the oxygen charge state is at least Q=4 in 
the region of the plasma torus, and therefore, in order to determine an 
upper limit to the energy loss of the energetic ions in the plasma of the 
inner magnetosphere, it will be assumed that the ions have reached their 
equilibrium charge state, Q=Z*. Again ions with 1 UeV/nuc are 
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Hgure 4.4 

a) The calculated lower limit of the oxygen ionization state as a function 
of radial distance from Jupiter. The values plotted represent the average 
distance at which a given ionization state is reached. The approximation 
was made in the calculation that the only process affecting the oxygen 
charge state is electron-impact ionization due to free electrons in the 
plasma of the inner magnetosphere. It was assiimed that the oxygen ions 
entering the region are singly iorazed. 

b) The oxygen g3rrofrequency for ions with the ionization state shown in 
panel a. A dipolar magnetic field was assumed. 
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considered (the energy loss is less for higher energies), so that QmO 
(4.10). The equation for the energy loss in a plasma is 


~ = -2.4x10 
dx 



i/s 


+ 9.13 


(4.11) 


(see e.g. Gloeckler at of. [1875]) where E is in MeV/nuc, x is the path 
length (cm). n« is the electron density (cm~^). and is the electron 
temperature (^K). Two regions will be considered; one near B Rj where 
n^wlO* cmT* [Warwick «f aL, 1979] and T.^xlO® K [Scudder «f at., 1961], 
and one near 6 Rj (in the hot torus) where n««2xl0* cm“® and T,«10® K 
(same references). Using equation (4.11) with the above assumptions and 
msasurements gives a lower limit to the time for a 1 MeV/nuc oxygen ion 
to lose 0.1 MeV/nuc near 0 Rj of *«2 years compared with 10-TBtroa*^9 
days. In the hot torus, the lower limit is ~35 days compared with 
10 Trtpon«"*3 days. Therefore, it is clear that the energetic ions do not 
lose a significant amount of their energy to the plasma in the inner 
magnetosphere. 


It has been shown in this section that the observed losses of the 
energetic ions inside 12 Rj can not be explained by geometiic absorption 
by lo, absorption by dust, or energy loss in plasma. Hence, it is assumed 
that pitch-angle scattering into the loss cone is the dominant loss 
mechanism in this region. 


4.6 Eiwriy Flow Rates 

The ultimate sink of particles that are pitch^angle scattered into 
the loss cone is the Jovian atmosphere. The total power that the 
energetic oxygen and sulfur ions deliver to the atmosphere will now be 
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calculated and coicpared with the power required to excite the aurora. 

The energy flow rate, or power, P, di.h^ered by energetic oxygen 
ions with magnetic moment greater than H across radius r is given by 

P>ll — p>||*2Tir*2h'VD‘I^ W (4.12) 

where p>y is the energy density of the oxygen ions (joules/ cn?), and the 
other quantities are defined following equation (4.6). In order to 
calculate P>y at as low a magnetic moment as possible and sUU be 
outside the region in which most particles arc lost, we have determined 
P>g at lO.i Rj. The quantity p>g was calculated from the 10.1 Rj 
spectrum in Figure 3.8, using the dash'dotted portion of the spectrum for 
magnetic moments less than 4.0x10^ MeV/nuc-G. As discussed in section 
3.5.3, the shape of this extrapolated curve was derived from the energy 
per charge spectrum, at the same M, measured by lilCP at '^IS Rj 
[Hamilton ef of.. 1861]. 

Figure 4.5 shows the values of P>n(10.lHj) calculated with equation 

(4.12) for several magnetic moment thresholds. The open symbols 

represent measurements based on the extrapolated portion of the 

spectrum. The uncertainties included in the error bars are those due to 

1 

uncertainties in abundance, event nunsber, And D. The uncertainty 

in the value of D at 10.1 Rj for different magnetic moments was estimated 
to be 50% (see Figure 4.2). An additional factor of 3 absolute uncertainty 
in D applies to all measurements. The power due to sulfur ions is 
represented by the X s 3 rmbols in Figure 4.5, and was determined from the 
oxygen P>n values using the S/0 ratios in Figure 3.5 and the relationship 
between the oxygen and sulfur spectral shapes (7s=To-f0.2) implied by 
that figure. The sulfur power is approximately equal to the oxygen power 
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ngure 4.5 

The inward energy flow rate across 10 Rj due to oxygen and suUixt ions 
with magnetic moment? greater than M; calculated using equation (4.12). 
The oxygen energy density at 10 Rj used in the calculation was 
determined from the 10.1 spectrum in Figure 3.8. The points with 
open circles were obtained using the extrapolated portion of the 
spectrum derived from LECP measurements. The energy flow rate due to 
sulfur ions was determined from the oxygen values using the sulfur to 
oxygen ratio in Figure 3.5 and the relationship between the oxygen and 
sulfur spectral shapes (7s = 7o40.2) implied by that flgiire. The energy 
flow rate across 10 Rj is approximately the power delivered to the Jovian 
atmosphere by ions lost between 10 and 8 Rj. The dotted curve 
represents the power delivered to the atmosphere by all oxygen and 
sulfui' ions lost from the magnetosphere inside 12 Rj. Also shown is the 
power required to produce the observed ultraviolet auroral emissions 
[Broadfoot cf ai., 1981]. 
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due to the fact that the sulfur ions are approximately half as abundant 
as oxygen ions, but have twice as much total energy (twice as many 
nucleons) at a given magnetic moment. 

The inward energy flow across 10 Rj due to >70 MeV/nuc>G oxygen 
and sulfur ions is seen to be ~1.5xl0** W. An interesting number to 
compare this with is the total outward energy flow due to the plasma 
ions. The number of plasma oxygen and sulfur ions diilusing outward 
across 10 Rj is «^5xl0^ ions/s, and their average total energy is 
approximately the corotation energy of an ion with mass 21 amu (the 
mean atomic mass of the dissociation products of SO 2 ) which is 1.7 keV. 
The outward energy flow is therefore ~1.4xl0*® W. Hence, the inward flow 
of >70 MeV/nuc>G oxygen and sulfur ions is comparable to the outward 
energy flow of the plasma ions. The energetic oxygen and sulfur inward 
energy flow increases toward lower magnetic moments, so that, including 
the lower energy ions and the energetic protons and electrons that also 
diffuse inward in this region, there is a substantial net energy flow into 
the inner magnetosphere. 

Since roost of the ions crossing 10 Rj are lost before reaching 8 Rj 
(Figure 4.3), the energy flow rate shown in the Figure 4.5 is. to flrst 
order, the power delivered to the atmosphere by ions lost between 10 and 
Rj. However, ions gain energy as they diffuse inward, so that the 
power from inside 8 Rj is not necessarily negligible. An estimate of the 
ratio of the energy flow crossing 6 Rj to that crossing 10 Rj is given by 

P>ll(BRl) „ 

P>m( 10R|) I>y(10Bj) IdORi) ‘ 

where l:(r)*‘U B(r) is the everege energy of ions with magnetic moment 
greater than U. Based on the inward flow rates of >4.0x10^ UeV/nuc-G 
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oxysen ions in Figure 4.3 and the measured magnetic Held strengths at 8 
and 10 Rj. P>h(8R;j)/P>u(10Rj)«0.S. The ion flow rate drops off steeply 
Inside 8 Rj. implying that most of the power, P>n(8Rj). is delivered to the 
atmosphere by ions lost between ~7.5 and 8 Rj, and that a relatively 
small amount is delivered from inside Rj. Similarly for 12 based 
on the inward flow rate of >1.7x10^ MeV/nuc-G ions in Figxire 4.3, 
P>h(12Ri)/P>ii(10Rj)mi.5. Outside 12 Rj losses are small. Hence, it is 
estimated that the total power delivered by sulfur and oxygen ions inside 
~12 Rj to the atmosphere is a factor of 6 = 2x(14*l.S*f0.5) larger than the 
oxygen energy flow rate across 10 Rj. This power is shown as the dotted 
curve in Figure 4.5. 

As discussed in Chapter 1, the most energetic phenomenon related 
to the inner magnetosphere is the aurora. The extreme ultraviolet 
auroral emission observed by the Voyager ultraviolet spectrometer 
requires a continuous power input of 1.2x10^^ W into the auroral zone 
[Broadfoot at., 1961]. This number is likely a lower limit because it 
was calculated assuming that all emission occurred high enough in the 
atmosphere to neglect extinction, whereas there is evidence that some of 
the emission originates deeper in the atmosphere where extinction is 
non-negligible [D. Shemansky. private communication]. The 
measurements shown in Figure 4.5 indicate that, if the required power is 
not substantially greater than 1.2x10^^ W, oxygen and sulfur ions with 
energies of «»100 keV/nuc are contributing a reasonable fraction of it. 
The figure also shows that lower energy ions are likely depositing 
additional power. As seen above, the bulk of the power delivered by the 
energetic oxygen and sulfiir ions comes from ions lost between ~7.5 and 
*^12.5 Rj. which, using the Voyager magnetic fleld models of Connemey ef 
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al. [1981], maps into a latitudinal 2 one in the atmosphere of 
This is consistent with the observed auroral zone width of roughly S” (see 
Chapter 1). and is within the range of latitudes at which aiu*orae are seen 
[Broadfoot et al., 1961]. 

As a final point, we note that the atmospheric penetration depths 
of heavy ions and electrons of equal energies are quite different as 
illustrated in Figure 4.6. Therefore, further analysis of the depth 
structure of the optical and ultiaviolet auroral emission may help 
diKtingiiiiih between excitation caused by electrons lost from the torus, 
and that caused by energetic oxygen and sulfur ions lost between 7.5 and 


12.5 Rj. 
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Rgure 4^ 


The column depth of penetration of electrons, protons, oxygen ions, and 
sulfur ions in Hg (the most abundant species in the upper Jovian 
atmosphere), as a fimcUou of their incident energies, assiiming the 
particles are incident normal to the atmosphere. The electron curve 
represents the extrapolated range given by 


R = 


^ 12.412 

Z 27 




[Koral and Cohen, 1965]. The proton, oxygen, and sulfur curves were 
calciUated using the range-energy tables of Northcliffe and Schilling 
(1970]. 


Column Depth of Penetrotion (cm*^) 
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Chapter 5 
Conclusions 

The principal conclusions ol this study of Z>2 ions with energies in 
the range 1 to 20 MeV/nuc between 5 and 20 Rj in the Jovian 
magnetosphere are as follows: 

1) The elenvental abundances of carbon, sodium, and sulfur relative to 
oxygen vary as a function of particle magnetic moment. H (*^p^VB), but 
are to first order the same throughout the magnetosphere at a given M. 
For 10® MeV/nuc-G ions (E»4 MeV/nuc at 10 Rj) the following abundance 
ratios are found: S/0=0-5, C/0=0.1, and Na/0«0.05. The sulfur to 
oxygen and sodium to oxygen ratios increase toward smaller M lmpl 3 nng 
that the sulfur and sodium energy spectra are softer than that of 
oxygen, whereas the carbon to oxygen ratio decreases toward smaller M 
Implying a hau’der carbon spectrum than that of oxygen. 

2) Between 6 and 17 Rj, there is a positive radial gradient in the phase 
'-pace density of energetip oxygen ions with magnetic moments in the 
range 2x10® to 2x10® MeV/nuc-G. The diffusive flow of the energetic ions 
is therefore inward in this region. 

3) The diffusion coefllcient upper limit at 9 Rj is ~10"® s“'. This limit, 
combined with the analysis of Voyager plasma observations by S^scoe et 
al. [1981], specifies an upper limit to the mass loading rate near lo of 
~10®® ions/s. 

4) A substantial fraction of the energetic oxygen ions are lost as they 
diffuse from 17 to 6 Rj, with the largest total niimber of particles lost 
between 8 and 12 Rj. The particle lifetime throughout the region is 
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within an order of magnitude of the strong pitch>angle diffusion lifetime. 
The location of the maxinium losses (6~12 Rj) indicates that geometric 
absorption by lo (5.9 Fcj) is not the primary loss mechanism for energetic 
oxygen ions. Also, order of magnitude estimates show that absorption by 
dust grains and energy loss in the plasma of the inner magnetosphere 
can not produce the observed losses. The loss mechanism is therefore 
postulated to be pitch-angle scattering into the loss cone. 

5) The total number of >70 MeV/nuc-G oxygen and sulfur ions diffusing 
inward across 10 Rj is ~10^ ion/s. This is a small number conqiared with 
lo’s source strength of 10*’— 10*® ions/s. The energy flow inward across 
10 Rj due to the same energetic ions, however, is '>-10^* W, which is 
approximately equal to the total outward energy flow at 10 Rj due to 
plasma ions. The energetic oxygen and sulfur inward energy flow 
increases toward lower magnetic moments, so that, including the lower 
energy ions and the energetic protons and electrons that also diffuse 
inward, there is a substantial net energy flow into the inner 
magnetosphere. 

6) The energetic oxygen and sulfur ions that are scattered into the loss 
cone during the diffusion process deposit their energy in the Jovian 
atmosphere. Based on the observed losses between 6 and 17 Rj, it is 
estimated that oxygen and sulfvir ions with magnetic moments greater 
than 70 MeV/nuc-G (E20.3 MeV/nuc at 10 Rj) deliver "-5x10** W to the 
Jovian atmosphere; a power comparable to the ~1.2xl0*® W required to 
produce the observed ultraviolet auroral emission. Additional power is 
hkeJb^ contributed by oxygen and sulfur ions with magnetic moments less 
than 70 MeV/nuc-G. Since most of the particles are lost between 7.5 and 
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12.5 Rj. the latitudinal zone in which the most energy is deposited in the 
atmosphere is ~67® to ~71®. 
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Appendix A 
Oeflnition of Tenm 

Brief definitions of several magneiospheric terms used in this 

thesis are given below. More details can be found in Roederer [1970] and 

in Schulz and Lanzerotti [1974]. 

bow shock — the shock-wave surface on the simward side of the 
magnetopause where the solar wind changes from a supersonic to a 
subsonic plasma. 

loss cone — the cone formed at any point along a magnetic field line by 
rotating a line at angle with respect to the field line about the 
field line. The angle O] is defined such that all particles with pitch 
angles less than aj will propagate far enough along the field line to 
encounter Jupiter's atmosphere before mirroring back toward the 
magnetic equatorial plane. The particles that encounter the 
atmosphere are absorbed and therefore lost from the 
magnetosphere. At 10 Rj. the loss cone at the magnetic equatorial 
plane has 01 ^ 1 °. 

magnetopause - the surface that bounds the magnetosphere, 

separating the domain of the inteiplanetry magnetic Held and 
plasma from that of the Jovian magnetic field. 

magnetosphere — the region of space in which Jupiter’s magnetic field is 
dominant 

MeUwoin L parameter (or simply L parameter) - a coordinate often used 
in magnetospheric physics to order observations of trapped 
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particles. Defined for a dipolar magnetic field as the distance from 
the center of the dipole that a particle's guiding center field line 
crosses the magnetic equatorial plane. L is defined in units of Dj, 
and is therefore dimensionless. In this thesis, L values were 
assigned only to observations made at or near the magnetic 
equatorial plane, and were therefore defined as the spacecraft - 
dipole center distance. (Jupiter's dipole center is offset ~0.1 Hj 
from the planet's center of mass). 

phase space density ~ the particle density in coordinate - momentum 
space (six dimensions). 

pitch angle — the angle between a particle's momentum vector and the 
magnetic field direction (i.e., pitch angles of 0° or 160** imply field- 
aligned flow). 

zenocentric - with respect to Jupiter’s center of mass, zeno- comes 
from Zeus. A table of the commonly used prefix and genitive forms 
of the planetary names is in Allen [1973, p.l39]. 
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Appendix B 

High Counting Rate Cahbrationa 

The high radiation levels of Jupiter's inner magnetosphere cause 
several instrumental problems in the LET system for which corrections 
must be made before absolute particle fluxes and energies can be 
determined. These include discriminator threshold changes due to 
baseline shifts at the output of the detector amplifiers, discriminator 
deadtime effects, and pulse pileup effects. In order to correct for the 
deadtinoe and threshold effects, a spare CRS instrument was calibrated 
using light pulses from randomly triggered light emitting diodes to 
simulate the ionization energy loss of charged particles in the detectors. 
The calibration and its results are described in detail elsewhere [Gehrels, 
1981]. This appendix summarizes the calibration results that are 
relevant to the present study and describes the pulse pileup effect. 


B.1 Discriminator Deadtimes 


Finite retrigger times for discriminator circuits cause a dead time 
in the instrument response that results in observed rates that are lower 
than the true rates. For the CRS instrument, the effect becomes 
noticeable above rates of ~10^ s~^ and significant above rates of 10^ s~^. 

The standard livetime correction formula for rate R is 


D _ n>i»d 

^ ” (I - 


(B.l) 


where r is the retrigger time of the discriminators, often assumed to be a 
constant. During the calibration it was determined, however, that, for 
the Z>2 rate and the detector smiles rates, t is a function of both the 
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rate and pulse height of the particles exceeding the discriminator 
threshold. It was therefore necessary to measure r. or equivalently 
Robaerrad ^ function of Rtni«. for all rates and conditions of interest. 

In Table B.l. measurements of Z>2oba a function of Z>2tnw are 
listed for 15» 30, and 90 HeV pulses in LI. In Figure 2.5 it is seen that 15 
MeV pulses are typical of those exceeding the Z>2 discriminator 
threshold outside of lo's orbit, and 30 and 90 MeV correspond to oxygen 
and sulfur events respectively, inside Io‘s orbit Although not shown in 
the figure, the LI pulse-height histogram at '^8 Rj has a single maximum 
at 30 MeV. The Z>2 livetime correction listed in Table 3.2 is the ratio 
Z>2tni«/ Z> 24 bt lor the observed LI pulse-height distribution in each 
region. 

The livetime correction formula for a coincidence rate such as LZ3 
is approximately given by 

° (1 - t.RuKI - t,Ru)(1°- - T2«Ri«) 

where Ri and ti are the singles rate and retrigger time of detector Li, and 
Rz >2 and T 2>2 are the Z>2 discriminator rate and retrigger time. 
Equation (B.2) is exact only in high fiux environments where all of these 
rates are much larger than the observed coincidence rate. During the 
calibration, the livetime correction factor, \ = l/Cl-TjRi), for each 
discriminator in the LZ3 coincidence equation (LLL2 L3 command state; 
14 off) was measured for the case where the particles of interest are 
oxygen and sulfur nuclei (large pulses) and the predominant events are 
protons (small pulses). 

Figiire B.1 shows fits to the measured Ai for 30 MeV (oxygen) and 
90 MeV (sulfur) pulses as a function of the rate. Rj. of MeV (proton) 
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Table B.l 


Observed vs True Z>2 Rate 
for Several LI Pulse Heights ^ 


Z>2 


true 


(10^ s"^) 


El-15 MeV 
(10^ s“^) 


Z>2 


obs 


El-30 MeV 
(10^ s‘^) 


El-90 MeV 

(10^ s"b 


5 

4.7 

4.7 

4.7 

7 

6.5 

6.5 

6.3 

10 

8.9 

8.8±0.2^ 

8.6 

15 

12.6 

12.3 

11.6 

20 

15.9 

15.5 

14.0 

25 

18.8 

18.3 

15.8 

30 

20.3 

21.0±1.0 

17.5 


E2-0, E3-0. 

^ Typical standard deviation of the measured Z>2 u,. 
in LET A. B. C. and D, 
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ngure R1 

Fits to the measiired LI livetimc correction factor, Aj, for the LZ3 rate 
for 30 Mev (oxygen) and 90 MeV (siilfur) pulses as a function of the rate. 
Ri, of MeV (proton) pulses. The error bars at the end of the ciirves 
represent the typical spread in the data about the fit Also shown are 
lines of constant of constant Tj in the equation Ai = 1/(1 -TjRj). 
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pulses. The error bars indicate the typical spresd in the data about the 
fit, and the dashed lines are lines of constant T|. The oxygen and sullui' 
curves are seen to be the same to within the ineasureinent uncertainty, 
so that, for a given Ri, the same LI correction factor will be applied to 
both event types. The results deviate significantly from lines of constant 
T, indicating the importance of directly roeasiiring tne correction factors. 

The measured As and for the above case were found to be 
consistent with Ag = Aa= 1.0 to within ~2% for all L2 and L3 rates. This is 
due primarily to the fact that the LETs use a strobed coincidence logic 
triggered by Ll (see e^. Gehrels [1061]). An event that exceeds the Li 
threshold and produces a strobe is counted in the LZ3 rate if it either 
triggers the L2 and L3 discriminators itself or comes at a time when 
these discriminators are up from a previous event. Hence, the L2 and L3 
discriminators cause very little deadtime in the LZ3 rate. It was also 
determined that Xz>sM 1.0, due both to the above effect and to the fact 
that approximately half of the events that exceed the Z>2 discriminator 
also cause an Ll singles rate count and are therefore already included in 
the correction. With these measurements and considerations, equation 
(B.2) simplifies to 


LZ3tpy* — Xx(Rj) LZ3ob > cr»«d • (B.3) 

The values of X} in the regions used in the analysis are listed in Table 3.2 
(the LZ3 livetime correction). 

Bw2 Diaciliianator Threshold Shifts 

fhgh counting rates cause negative shifts in the baseline voltage 
level at the output of the detector anpliflers. Since discriminator 
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tiirssholds are let for a ipecifled level above the nominal baseline, the 
baseline shift produces an increase in the effective threshold, so that, for 
instance, at high counting rates it takes a larger energy deposit to 
exceed a given threshold than at low counting rates. During the 
calibration, the shift in the LI Z>2 threshold was measured as a function 
of the LI pulse height and Z>2 rate. The effective carbon, oxygen, and 
sulfur Z>2 thresholds, including energy loss in the window, were 

then calculated from the measured thresholds and are listed in Table B.2. 
The effective oxygen thresholds used in the analysis are listed in Table 
3.2 


B.3 Pulse pileiqt 

As fluxes increase, the probability that two particles are coincident 
in a detector within the instrumental resolving time increases, degrading 
the pulse-height data from that detector. For large oxygen pulses in the 
detector, the effect of a high rate of coincident small proton pulses is to 
broaden the oxygen pulse*height distribution. Oxygen events in 
coincidence with other large pulses, however, will be displaced conq^lelely 
cut of the normal oxygen distribution. The fraction displaced is 
Tp>|» Z>2tnw, where is the effective resolving time for pulse>height 
analysis (~20 a«s). In all regions, particles stopping in L3 were analyzed 
using only their L2 and L3 pulse-height information. The maximum 
counting rate of large oxygen and sulfur pulses in these two detectors is 
on the order of 10^ s"*^ (L2 near lo), so that pulse pileup affects <2% of 
the Z>2 particles stopping in L3. In the 4.9, 5.3, and 5.9 I<j regions 
however, particles stopping in L2 were used to calculate the lower energy 
portions of the spectra, and for these events the LI p\ilse height is 
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Table B.2 


Effective Z>2 Discriminator Thresholds 


c ^ 

Pulse *'t-eff 

- , Height Carbon Oxygen Sulfur 

(10'^ s"M (MeV) (MeV/nuc) (MeV/nuc) (MeV/nuc) 


6.6 

7.5 

1.2 

1.1 

0.9 

11.6 

7.5 

1.3 

1.1 

0.9 

19.7 

7.5 

1.3 

1,1 

0.9 

6.5 

15.4 

1.3 

1.1 

0.9 

11.6 

15.4 

1.4 

1.2 

0.9 

19.7 

15.4 

1.5 

1.3 

1,0 

6.6 

30.8 

1.3 

1.2 

0.9 

11.6 

30,8 

1.5 

1.3 

1.0 

19.7 

30.8 

1.7 

1,4 

1.1 

6.6 

92.3 

1.4 

1.2 

1,0 

11.6 

92.3 

1.7 

1.4 

1.1 


t 


The nominal thresholds for carbon, oxygen, and sulfur 
are 1.2, l.I. and 0.9 MeV/nuc respectively. 
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required lor elemental identillcation. The counUii^. rate of large pulees in 
Ll ie Z>2tra» wfdch la ~i.SxlO^ a*** in theae regiona, ao that piilae pileup 
affecta approxiznaielj 30X of the Z>2 eyenta. A more preciae eatixnate of 
the amount of pulae pileup in each region can be determined from the 
data by cong>aring the energy apectium of partinlea atopping in L3 
obtained uaing Ll in the elemental identification proceM with one in 
which only L2 and L3 were uaed. The ratio between the reairlting apectra 
ia energy independent over the three^detector energy range and indicates 
the fraction of the events in Ll displaced out of the oxygen distribution. 
The ratio was therefore used as the pulse pileup correction factor tor the 
two*detector part of the spectrum 
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